arXiv:1502.01586v2 [astro-ph.IM] 17 Jul 2015 


Astronomy & Astrophysics manuscript no. HFIDPC © ESO 2015 

July 20, 2015 


Planck 2015 results. VII. High Frequency Instrument data 
processing: Time-ordered information and beams 

Planck Collaboration: R. Adam’*, P. A. R. Ade*’, N. Aghanim^^, M. Amaud®, M. Ashdown®-^, J. AumonP^, C. Baccigalupi*', A. J. Banday^®-®, 
R. B. Barreiro®, N. Bartolo’*’^', E. Battaner®'-^’, K. Benabed^®’*®, A. BenoiP^, A. Benoit-Levy’’’^® ®, J.-P. Bernard®®’^, M. Bersanelli^*’"*^, 

B. BertincourP^, P. Bielewicz’**’*', J. J. Bock® *', L. Bonavera®, J. R. Bond*, J. Borrill**’*^, F. R. BoucheP® *'*, F. Boulanger^^, M. Bucher*, 

C. Burigana'*^-’'*-'**, E. Calabrese*’, J.-F. Cardoso’'*•*■*^ A. Catalano’*®, A. Challinor^*® *’, A. Chamballu®-*^-**, R.-R. Chary®, H. C. Chiang’’.^, 

P. R. Christensen’**®, D. L. Clements®, S. ColombE*****, L. P. L. Colombo’*-®, C. Combet’*, F. Couchot®, A. Coulais®, B. P. Crill®-** *, 

A. Curto®-*’-®, F. Cuttaia'*"', L. Danese**, R. D. Davies®*, R. J. Davis®*, P. de Bemardis***, A. de Rosa'*"*, G. de Zotti***-**, J. Delabrouille*, 
J.-M. Delouis*®-*®, F.-X. Desert®, J. M. Diego®®, H. Dole®*-®, S. Donzelli®, O. Dore®’-**, M. Douspis*®, A. Ducout®®-®*, X. Dupac*®, 

G. Efstathiou®*, F. Eisner”-®®-*®, T. A. EnBlin’®, H. K. Eriksen®®, E. Falgarone®’, J. Fergusson*’, F. Finelli®-'*®, O. Forni®®-®, M. Frailis®*, 

A. A. Fraisse’®, E. Franceschi®, A. Frejsel’®, S. Galeotta®*, S. Galli®'*, K. Ganga*, T. Ghosh®®, M. Giard®®-®, Y. Giraud-Heraud*, E. Gjerlpw®®, 

J. Gonzalez-Nuevo*’-®®, K. M. Gorski®’-®*, S. Gratton®®-®*, A. Gruppuso®, J. E. Gudmundsson’®, F. K. Hansen®®, D. Hanson’®-®’-*, 

D. L. Harrison®*-®®, S. Henrot-Versille®®, D. Herranz®®, S. R. Hildebrandt®’-**, E. Hivon®®-*®, M. Hobson®, W. A. Holmes®’, A. Homstrup*®, 

W. Hovest’®, K. M. Huffenberger”, G. Hurier®®, A. H. Jaffe®*, T. R. Jaffe®®-®, W. C. Jones’®, M. Juvela’”*, E. Keihanen’®, R. Keskitalo**, 

T. S. Kisner’*, R. Kneissl*"*-’, J. Knoche’®, M. Kunz*®-®®-’, H. Kurki-Suonio’®-'*®, G. Lagache®-®®, J.-M. Lamarre®’, A. Lasenby®-®®, M. Lattanzi’®, 
C. R. Lawrence®’, M. Le Jeune*, J. P. Leahy®*, E. Lellouch®*, R. Leonard!*®, J. Lesgourgues®’ **, F. Levrier®’, M. Liguori’*-®*, P. B. Lilje®®, 

M. Linden-Vprnle*®, M. Lopez-Caniego*®-®®, P. M. Lubin’®, J. F. Maci'as-Perez’*, G. Maggio®*, D. Maino**-®®, N. Mandolesi®**-’®, A. Mangilli®®-®®, 
M. Maris®*, P. G. Martin*, E. Martmez-Gonzalez®®, S. Masi*®, S. Matarrese’*-®*-**, P. McGehee®’, A. Melchiorri*®-®’, L. Mendes*®, 

A. Mennella**-®®, M. Migliaccio®*-®®, S. Mitra®®-®’, M.-A. Miville-Deschenes®®-*, A. Moneti®®, L. Montier®®-®, R. Moreno®*, G. Morgante®®, 

D. Mortlock®*, A. Moss**, S. Mottet®®, D. Munshi*’, J. A. Murphy”, P. Naselsky’®-**, F. Nati’®, P. Natoli’®-*-®®, C. B. Netterfield**, 

H. U. Nprgaard-Nielsen*®, F. Noviello®*, D. Novikov’®, 1. Novikov’®-’®, C. A. Oxborrow*®, F. Paci**, L. Pagano*®-®’, F. Pajot®®, D. Paoletti®®-®®, 

F. Pasian®*, G. Patanchon*, T. J. Pearson**-®’, O. Perdereau®®, L. Perotto’*, F. Perrotta**, V. Pettorino*®, F. Piacentini*®, M. Piat*, E. Pierpaoli’*, 

D. Pietrobon®’, S. Plaszczynski®®, E. Pointecouteau®®-®, G. Polenta*-®’, G. W. Pratt®®, G. Prezeau**-®’, S. Prunet®®-*®, J.-L. Puget®®, J. P. Rachen*®-’®, 
M. Reinecke’®, M. Remazeilles®*-®®-*, C. Renault’*, A. Renzi*’-®*, 1. Ristorcelli®®-®, G. Rocha®’-**, C. Rosset*, M. Rossetti**-®®, G. Roudier*-®’-®’, 
M. Rowan-Robinson®*, B. Rusholme®’, M. Sandri®®, D. Santos’*, A. Sauve®®-®, M. Savelainen’®-®®, G. Savini*®, D. Scott’®, M. D. Seiffert®’-**, 

E. P. S. Shellard*’, L. D. Spencer*’, V. Stolyarov®-®®-*®, R. Stompor*, R. Sudiwala*’, D. Sutton®*-®®, A.-S. Suur-Uski’®-®®, J.-F. Sygnet®®, 

J. A. Tauber*®, L. Terenzi*’-®®, L. Toffolatti*’-®®-®®, M. Tomasi**-®®, M. Tristram®®, M. Tucci*®, J. Tuovinen*®, L. Valenziano®®, J. Valiviita’®-®®, 

B. Van Tent”, L. Vibert®®, P. Vielva®®, F. Villa®®, L. A. Wade®’, B. D. Wandelt®®-*®-”, R. Watson®*, 1. K. Wehus®’, D. Yvon*®, A. Zacchei®*, 

A. Zonca’® 

(Affiliations can be found after the references) 

Submitted to A&A February 8, 2015; Accepted July 15, 2015 


Abstract 

The Planck High Frequency Instrument (HFI) has observed the full sky at six frequencies (100, 143, 217, 353, 545, and 857 GHz) in intensity and 
at four frequencies in linear polarization (100, 143, 217, and 353 GHz). In order to obtain sky maps, the time-ordered information (TOI) containing 
the detector and pointing samples must be processed and the angular response must be assessed. The full mission TOI is included in the Planck 
2015 release. This paper describes the HFI TOI and beam processing for the 2015 release. HFI calibration and map making are described in a 
companion paper. The main pipeline has been modified since the last release (2013 nominal mission in intensity only), by including a correction 
for the nonlinearity of the warm readout and by improving the model of the bolometer time response. The beam processing is an essential tool 
that derives the angular response used in all the Planck science papers and we report an improvement in the effective beam window function 
uncertainty of more than a factor of 10 relative to the 2013 release. Noise correlations introduced by pipeline filtering function are assessed using 
dedicated simulations. Angular cross-power spectra using data sets that are decorrelated in time are immune to the main systematic effects. 
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1. Introduction: a summary of the HFI pipeline 

This paper, one of a set associated with the 2015 Planck^ data 
release, is the first of two that describe the processing of the 
data from the High Frequency Instrument (HFI). The HFI is one 
of the two instruments on board Planck, the European Space 
Agency’s mission dedicated to precision measurements of the 
cosmic microwave background (CMB). The HFI uses cold op¬ 
tics (at 4K, 1.6 K, and 100 mK), filters, and 52 bolometers 
cooled to 100 mK. Coupled to the Planck telescope, it enables 
us to map the continuum emission of the sky in intensity and 
polarization at frequencies of 100, 143, 217, and 353 GHz, and 
in intensity at 545 and 857 GHz. Paper 1 (this paper) describes 
the processing of the data at the time-ordered level and the mea¬ 
surement of the beam. Paper 2 (Planck Collaboration VIII 2015) 
describes the HFI photometric calibration and map making. 

The HFI data processing for this release is very similar to 
that used for the 2013 release (Planck Collaboration VI 2014). 
Figure 1 provides a summary of the main steps used in the pro¬ 
cessing, from raw data to frequency maps both in temperature 
and polarization. 

First, the telemetry data are converted to TOI (time ordered 
information). The TOI consists of voltage measurements sam¬ 
pled at 180.3737 Hz for each of the 52 bolometers, two dark 
bolometers, 16 thermometers, and two devices (a resistance and 
a capacitance) that comprise the HFI detector set. The TOI is 
then corrected to account for nonlinearity in the analogue-to- 
digital conversion (see Sect. 2). Glitches (cosmic ray impacts on 
the bolometers) are then detected, their immediate effects (data 
around the maximum) are flagged, and their tails are subtracted. 
A baseline is computed in order to demodulate the AC-biased 
TOI. A second-order polynomial correction is applied to the de¬ 
modulated TOI to linearize the bolometer response. The minute- 
scale temperature fluctuations of the 100 mK stage are subtracted 
from the TOI using a combination of the TOI from the two dark 
bolometers. Sharp lines in the temporal power spectrum of the 
TOI from the influence of the helium Joule-Thomson ("^He-JT) 
cooler (hereafter called 4-K lines) are removed with interpola¬ 
tion in the Fourier domain (see Sect. 3.3). The finite bolome¬ 
ter time responses are deconvolved from the TOI, also in the 
Fourier domain. For this release, the time response consists of 
four to seven thermal time constants for each bolometer. Several 
criteria based on statistical properties of the noise are used to re¬ 
ject the stable pointing periods (hereafter called rings) that are 
non-stationary (see Table 1). A subtractive jump correction is 
applied; it typically affects less than 1 % of the rings and the am¬ 
plitude of the jumps exceeds a tenth of the TOI rms in less than 
0.1 % of the rings. 

At this point, the TOI is cleaned but not yet calibrated. The 
beam is measured using a combination of planet observations 
for the main beam and GRASP physical optics calculations^ for 
the sidelobes (see Sect. 4). The focal plane geometry, or the rel¬ 
ative position of bolometers in the sky, is deduced from Mars 
observations. The 545 and 857 GHz channels are photometri¬ 
cally calibrated using the response to Uranus and Neptune. The 
lower frequency channels (100, 143, 217, and 353 GHz, called 

* Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from Prance and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 

^ TICRA, http://www.ticra.com/products/software/grasp 


the “CMB” channels) are calibrated with the orbital CMB dipole 
(i.e., the dipole induced by the motion of the Lagrange point L2 
around the Sun). 



1 


electronic 





& thermal 



Fourier transform 


responses 


beam 

1 


§3.4 


§-1 


HPR 

qualification 


computation 

of the data §5 



focal plane 
reconstruction 


computation of 
noise TOI 


Figure 1. Schematic of the HFI pipeline. The left part of the 
schematic involves TOI and beams (this paper), while the upper- 
right part represents the map making steps (Paper 2). Ancillary 
maps are composed of zodiacal light templates (ZL) and po¬ 
larization band-pass mismatch (BPM) maps. ADC = analogue- 
to-digital converter (see Sect. 2). HPR = HealPix ring (see 
Paper 2). “Beam products” refers to beam transfer functions, B(. 
Blue: changes in this release with section numbers correspond¬ 
ing to this paper. Yellow: released data products. 


In this paper, we describe the changes made to the processing 
since the 2011 and 2013 papers (Planck HFI Core Team 201 lb; 
Planck Collaboration VI 2014). Section 2 gives a view of a step 
that has been added to the beginning of the pipeline, namely the 
correction for ADC (analogue-to-digital converter) nonlineari¬ 
ties. This step proves to be very important for the quality of the 
CMB data, especially at low multipoles. Section 3 deals with the 
addition of more long time constants in the bolometer response 
and other TOI. Section 4 presents refined beam measurements 
and models. Some consistency checks are reported in Sect. 5. 
The public HFI data products are described in Appendix A. 

Paper 2 describes the new scheme for the CMB dipole cali¬ 
bration (the Planck orbital dipole calibration is used for the first 
time) and the sub millimetre calibration on planets. Paper 2 also 
describes the polarized map making, including the derivation of 
far sidelobes, and zodiacal maps, as well as polarization correc¬ 
tion maps due to bandpass mismatch: a generalized least-squares 
destriper is used to produce maps of the temperature and the 
two linear polarization Stokes parameters Q and U. About 3000 
maps are obtained by splitting the HFI data into different sub¬ 
sets by, e.g., time period or detector sets. Consistency checks are 
performed in order to assess the fidelity of the maps. 

2. ADC correction 

Planck Collaboration VI (2014) reported that the HFI raw data 
show apparent gain variations with time of up to 2% due to 
nonlinearities in the HFI readout chain. In the 2013 data re¬ 
lease (Planck Collaboration VIII 2014) a correction for this sys- 
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tematic error was applied as an apparent gain variation at the 
map making stage. The 2013 maps relied on an effective gain 
correction based on the consistency constraints from the recon¬ 
structed sky maps, which proved to be sufficient for the cosmo¬ 
logical analysis. 

For the 2015 data release we have implemented a direct ADC 
correction in the TOI. In this section, we describe the ADC effect 
and its correction, and its validation through end-to-end simula¬ 
tions (see also Sect. 5.4). Internal checks of residual ADC non¬ 
linearity are shown in Sect. 5. 

2.1. The ADC systematic error 

The HFI bolometer electronic readout (Planck Collaboration ES 
2015) includes a 16-bit ADC, the flight-qualified version of 
Maxwell Technologies model 7809LP having a very loose tol¬ 
erance on the differential nonlinearity (the maximum deviation 
from one least significant bit, LSB, between two consecutive lev¬ 
els, over the whole range), specified to be not worse than one 
LSB. The implications of this feature for HFI performance had 
not been anticipated, and it did not produce any detected effect 
in ground test data. However, it proved to be a major system¬ 
atic effect impacting the flight data. A wide dynamic range at 
the ADC input was needed to both measure the CMB and the 
sky foregrounds, and properly characterize and remove the tails 
of glitches from cosmic-rays. Operating HFI electronics with the 
necessary low gains increased the effects of the ADC scale errors 
on CMB data. 

We have developed a method that reduces the ADC effect 
on the angular power spectra by more than a factor of 10 for 
most bolometers. There are three main difficulties in making this 
correction. 

- The chip linearity defects were characterized with insuffi¬ 
cient precision before launch. We have designed and run a 
specific campaign to map the ADC nonlinearity from flight 
data (Sect. 2.2). 

- Each HFI data sample in the TOI is the sum of 40 consec¬ 
utive ADC fast samples , corresponding to half of the mod¬ 
ulation cycle. The full bandwidth of the digital signal is not 
transmitted to the ground. The effective correction of TOI 
samples due to ADC defects requires the knowledge of the 
shape of the fully-sampled raw data at the ADC input. A 
shape model is built from the subset of fully-sampled down¬ 
loaded data, transmitted to the ground at a low rate (80 suc¬ 
cessive fast samples every 101.4 s for every detector), here¬ 
after called “the fully-sampled raw data.” 

- The “^He-JT 4-K cooler lines in the TOI (Sect. 3.3) result 
from a complex parasitic coupling. Much of this parasitic 
signal is a sum of 20 Hz harmonics, synchronous with the 
readout clock, with nine slow sample periods fitting exactly 
into two parasitic periods. Capture of a sequence of 360 raw 
signal samples would have allowed a direct reconstruction 
of the full patterns of this parasitic signal. However, the short 
downloaded sequences of 80 samples always fall in the same 
4-K phase range, allowing us to cover only 2/9 of the full 
pattern. To properly model the signal at the ADC input, one 
must decipher this parasitic signal over its full phase range. 
The present model relies both on the full sampling subset 
and on the 4-K lines measured in the TOI (see Sect. 3.3 and 
references therein). 

The signal model, including the 4-K line parasitic part, is de¬ 
scribed in Sect. 2.3. 


2.2. Mapping the ADC defects 

The defects of an ADC chip are fully characterized by the in¬ 
put levels corresponding to the transitions between two consec¬ 
utive output values (known as digital output code, or DOC). An 
ADC defect mapping is usually run on a dedicated ground test 
bench. We made this measurement on two spare flight chips to 
understand the typical behaviour of the circuit in the range rele¬ 
vant to the flight data. These measurements revealed a 64-DOC 
periodic pattern, precisely followed by most of the DOC, ex¬ 
cept at the chunk boundaries, which allowed us to build a first, 
approximate defect model with a reduced number of parame¬ 
ters (Planck Collaboration VIII 2014). Such behaviour is under¬ 
stood from the circuit design, where the lowest bits come from 
the same components over the full ADC scale. However, because 
each ADC has a unique defect pattern, data from these ground 
tests could not be used directly to correct the flight data. 

The parameters for the on-orbit chips were extracted for each 
HFI bolometer using data samples of the thermal signal, effec¬ 
tively a Gaussian noise input. These ADC-dedicated flight data, 
herafter called “warm data,” were recorded during the 1.5 years 
of the LFI extended mission, between February 2012 and August 
2013. During this period the bolometer temperature was stable 
at about 4 K, a temperature at which the bolometer impedance is 
low, giving no input signal or parasitic pickup on the ADCs apart 
from a tunable offset and Gaussian noise with rms values around 
20 equivalent mean LSBs. The defect mapping was obtained by 
inverting the histograms of the accumulated fully-sampled raw 
data, as explained below. 

The warm data consist of large sets of ntot counts of 80 sam¬ 
ple raw periods taken in stable conditions with different bias cur¬ 
rents and/or compensation voltages on the input stage tuned to 
mostly sample the central area of the ADC. We denote the signal 
value at the ADC input v, and v, is the input level correspond¬ 
ing to the transition between DOC (i - 1) and DOC i. We call 
the probability density function of the k* sample of the 80 sam¬ 
ples p^(v). This is assumed to be a Gaussian distribution with a 
mean depending on k and on the data set, and a variance cr^ 
depending only on the data set. 

Every set k is a series of {n,j histogram bin statistics, follow¬ 
ing a multinomial distribution. The mean values, {«,), are linked 
to and (v,) by the following set of equations: 

X Vi+\ 

p’'(v)dv. (1) 

A maximum likelihood technique maps the {«,) to the {«,). To 
extract (v,) from Eq. (1) one needs an estimate of and cr^. The 
equations are solved by recurrence, starting from the most pop¬ 
ulated bin. Since the solution depends on and cr^, normally 
one would need to know their true value. In fact, this system has 
specific properties that allow recovery of the correct ADC scale 
without any other information. 

- Its solution (v,) is extremely sensitive to p*^. An incorrect in¬ 
put value gives unphysically diverging ADC step sizes; this 
allows us to choose the p^ that gives the most stable solu¬ 
tion. This prescription makes the (v,) solution in the popu¬ 
lated part of the histrogram independent from the hypothesis 
on p^. 

- The step size solution is proportional to cr, so the solutions 
for different data sets can be intercalibrated to fit a common 
ADC scale, which allows the unknown cr to be eliminated, 
provided all 80 samples of a data set have the same noise 
value. 
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Figure 2 shows a typical example of the maximum likeli¬ 
hood result, before applying the periodic feature constraints on 
the solution. This method is based on the strong assumption 
that the noise is Gaussian. Some data showing small tempera¬ 
ture drifts have not been used. With this method, the ADC step 
sizes are nearly independently measured, so there is a random- 
walk type of error on the DOC positions, which is limited by 
the 64 DOC periodicity constraint. The 64 DOC pattern is ob¬ 
tained from the weighted average of the values obtained in the 
+512 DOC around mid-scale, except for the first 64 DOC above 
mid-scale, which do not follow the same pattern. 

The main step is found at mid-scale. For all channels, there 
are more than 10^ samples per DOC histogram bin in the small 
+512 output code range around mid-scale that is explored by 
flight data. Outside this region, the smaller numbers of samples 
lead to bigger drifts in the likelihood solution. 

Both residual distributions and simulations give an estimate 
of the precision of the present defect recovery below 0.03 LSB 
for any DOC over this range. Systematic errors on large DOC 
distances, not taken into account in this model, are smaller than 
0.2 LSB over a range of 512 LSB. At this level, we see residuals 
due to violation of the rms noise (cr) stability assumption, and we 
are working on an improved version of the ADC defect recovery 
to be included in future data releases. 



32000 32200 32400 32600 32800 33000 33200 33400 
Output code 


Figure 2. An example of one ADC defect mapping around the 
mid-scale. The useful range +512 around mid-scale is shown 
with vertical lines. 


2.3. Input signal model 

The signal at the input of the ADCs is the sum of several compo¬ 
nents, including the modulation of the bolometer, noise from the 
bolometer, and the electronics, and bolometer voltage changes 
due to the astrophysical signal. The largest in amplitude is the 
modulation of the bolometer voltage with the combination of a 
triangle- and a square-wave signal, fed through bias capacitors. 
Most of the modulation is balanced (i.e., the voltage amplitude 
is reduced to nearly zero) by subtracting a compensating square- 
wave signal prior to digitization. 

Spurious components also exist, such as the 4-K parasitic 
signal and various electronic leakage effects. Both are well de¬ 
scribed by slowly varying periodic patterns, assumed to be stable 
on a one-hour time scale. 


The signal model used here is based on the linearity of the 
bolometer chain. It is a steady-state approximation of the signal 
shape produced by constant optical power on the bolometer. It is 
given by 

d(t) = PxG,,^(t) + 0(t), (2) 

where P is proportional to the sky signal, and the raw gain Graw 
and the offset O are periodic functions of time. The raw gain 
period is the same length as the readout period, and the offset 
period is equivalent to the 4-K cooler period. 

The parameters of the model have been extracted and 
checked over the whole mission from a clean subset of fully- 
sampled raw data (particle glitches and planet-crossings ex¬ 
cluded). For each bolometer, Graw is given by a set of 80 num¬ 
bers, assumed to be stable over the full mission. The offsets are 
given ring by ring as sets of 360 numbers. 

Figure 3 shows an example of the raw gain Graw, for bolome¬ 
ter 100-la. This function represents the fast time scale voltage 
across the bolometer. The detailed shape is due to the modulat¬ 
ing bias current and to well-understood filters acting on the sig¬ 
nal in the readout chain (for details, see Section 4.1 of Lamarre 
et al. (2010)). The signal model (Eq. 2) assumes that every pair 
of telemetered fast data samples is the integral of one half of this 
function, linearly scaled by the input sky signal. 

Figure 4 shows the evolution of the constant term over the 
mission, for the 143-6 bolometer, due to the spurious signals in 
the fully-sampled raw data window. 



Figures. Raw gain for the 100-la bolometer. The horizontal 
scale is fast sample index, while the vertical scale is arbitrary. 

Electronics leakage is readout-synchronous and is therefore 
taken into account by the present model. This is not the case for 
the 4-K cooler parasitic signal. The available information in the 
TOI does not allow reconstruction of the shape of the 4-K line 
at the ADC input without ambiguities. Eortunately, this parasitic 
signal is, in most cases, dominated by a few harmonics that we 
constrain from the combined analysis of the constant term in 
the fully-sampled raw data and the 4-K-folded harmonics in the 
TOI. A constrained method is used to extract amplitudes and 
phases ring per ring for the 20, 160, and 200 Hz components. 
Only these three harmonics are used in the present model. The 
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Figure 4. Constant term drift for the 143-6 bolometer, relative to 
ring 1000. The colour code goes from black to red over the full 
mission. 


effect of this approximation is checked by using both internal 
and global tests. Improved 4-K harmonic estimators will be in¬ 
cluded in future data releases. 

2.4. ADC correction functions 

The input shape model from Eq. 2 is used, along with the ADC 
DOC patterns and the 5 phase value, to compute, ring by ring, the 
error induced by the ADC as a function of the summed 40 sam¬ 
ples of the TOI. A set of input power values {P} is used to sim¬ 
ulate TOI data: first,with the input model that includes the 4-K 
line model and the ADC non linearity model; and, second, with 
a parasite-free input model and a linear ADC. We identify real 
TOI data with the first result and draw from these simulations 
the correction functions that are then applied to the real data. 

The combination of the two readout parities^ with different 
possible phases relative to the 4-K parasitic signal produces 18 
different correction functions that are determined for each data 
ring and applied to every TOI sample. Figure 5 shows an exam¬ 
ple. 

2.5. Simuiation of the ADC effect 

The measured ADC defects are included in the HFI end-to-end 
simulations (Sect. 5.4), run at the fully-sampled raw data level. 
This allows us to assess some consequences of the ADC non¬ 
linearity on the data. For instance, the gain variations are well- 
reproduced. Figure 6 shows an example of survey-difference 
angular power spectra. The systematic effect coming from the 
ADC behaves as l/{!^ at low and is well suppressed by the 
ADC correction. This is an important check on the consistency 
of the ADC processing. 


^ The science data samples alternate between positive and negative 
parity corresponding to the positive and negative part of the readout 
modulation cycle. 



ADU 

Figures. Example ADC correction functions for both parities 
(upper and lower curve sets) and two different phases (in black 
and red) relative to the '^He-JT cooler cycle. 



Figure 6. Survey-difference angular power spectra for the 100- 
4b bolometer. The 2015 release data (blue curve) are compared 
with simulated data containing ADC nonlinearities. Green and 
red curves show simulations produced, with and without the 
ADC correction, respectively. 


3. TOI processing 

After the ADC correction, the software modules used to process 
the 2015 TOI are identical to those used for 2013, up to the stage 
of clean TOIs ready for selection, calibration, and map making 
(see Fig. 1). However, some of the modules’ input parameters 
have been fine-tuned to better control residual systematic errors 
present in the 2013 data. We describe these tunings here, and 
later assess their accuracy in Sect. 5. 

3.1. Pointing and focai-piane reconstruction 

Satellite attitude reconstruction is the same for both Planck in¬ 
struments in the 2015 release and is described in the mission 
overview paper (Planck Collaboration I 2015). The major im¬ 
provement since 2013 is the use of solar distance and radiometer 
electronics box assembly (REBA) thermometry as pointing error 
templates that are fitted and corrected. 
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As in Planck Collaboration VI (2014), we determine the lo¬ 
cation of the HFI detectors relative to the satellite Foresight us¬ 
ing bright planet observations. Specihcally, we use the hrst ob¬ 
servation of Mars to dehne the nominal location of each detec¬ 
tor. These locations do not correspond exactly to the physical 
geometry of the focal plane, since they include any relative shift 
induced by imperfect deconvolution of the time-response of the 
detectors, as described in Sect. 3.4.2. After this initial geomet¬ 
rical calibration, the reconstruction is monitored by subsequent 
planet observations. We hnd that the detector positions are stable 
to < 10" with an rms measurement error of about 1". 

3.2. Cosmic ray deglitching 

The deglitching method, described in Planck Collaboration X 
(2014), consists of flagging the main part (with S/N > 3.3) of 
the response to each cosmic ray hit and subtracting a tail com¬ 
puted from a template for the remaining part. The flagged part is 
not used in the maps. The method and parameters are unchanged 
since 2013. As described by Planck Collaboration X (2014) and 
Catalano et al. (2014), three main populations of glitches have 
been identified: (i) short glitches (with a peaked amplitude dis¬ 
tribution), due to the direct impact of a cosmic particle on the 
grid or the thermometer; (ii) long glitches, the dominant popu¬ 
lation, due to the impact of a cosmic particle on the silicon die, 
or support structure of the bolometer’s absorber; (iii) and slow 
glitches, with a tail similar to the long ones and showing no fast 
part. The physical origin of this last population is not yet under¬ 
stood. 

The polarization-sensitive bolometers (PSBs) are paired, 
with two bolometers (called a and b) sharing the same hous¬ 
ing (Jones et al. 2003). The dies are thus superimposed and most 
of the long glitches seen in one detector are also seen in the other. 
A flag is computed from the sum of the a and b signal-subtracted 
TOI after each has been deglitched individually. This new flag is 
then included in the total flag used for both the a and b bolome¬ 
ters. 

For strong signals, the deglitcher threshold is auto-adjusted 
to cope with source noise, due to the small pointing drift dur¬ 
ing a ring. Thus, more glitches are left in data in the vicinity of 
bright sources (such as the Galactic centre) than elsewhere. To 
mitigate this effect near bright planets, we flag and interpolate 
over the signal at the planet location prior to the TOI process¬ 
ing. While a simple linear interpolation was applied in the first 
release (Planck Collaboration XII2014), an estimate of the back¬ 
ground signal based on the sky map is now used to replace these 
samples. For the 2015 release, this is done for Jupiter in all HFI 
frequency bands, for Saturn at v > 217GHz, and for Mars at 
V > 353 GHz. 

Nevertheless, for beam and calibration studies (see Sect. 4 
and Paper 2), the TOI of all planet crossings, including the 
planet signals, is needed at all frequencies. Hence, a separate 
data reduction is done in parallel for those pointing periods 
and bolometers. For this special production, the quality of the 
deglitching has been improved with respect to the 2013 data 
analysis (see Appendix B). 

3.3. ‘^He-JT cooler pickup and ring selection 

Planck scans a given ring on the sky for roughly 45 min before 
moving on to the next ring (Planck Collaboration I 2014). The 
data between these rings, taken while the spacecraft spin-axis is 
moving, are discarded as “unstable.” The data taken during the 


intervening “stable” periods are subjected to a number of statis¬ 
tical tests to decide whether they should be flagged as unusable 
(Planck Collaboration VI 2014). This procedure continues to be 
adopted for the present data release. Here we describe an addi¬ 
tional selection process introduced to mitigate the effect of the 
4-K lines on the data. 

The ^He-JT cooler is the only moving part on the spacecraft. 
It is driven at 40 Hz, synchronously with the HFI data acqui¬ 
sition. Electromagnetic and microphonic interference from the 
cooler reaches the readout boxes and wires in the warm service 
module of the spacecraft and appears in the HFI data as a set of 
very narrow lines at multiples of 10 Hz and at 17 Hz (Planck 
Collaboration VI 2014). The subtraction scheme used for the 
2013 release, used here as well, is based on measuring the 
Fourier coefficients of these lines and interpolating them for the 
rings that have one of the harmonics of their spin frequency very 
close to a line frequency - a so-called “resonant” ring. However, 
it was noticed that the amplitude of the lines increased in the last 
two surveys. Therefore, as a precautionary step, resonant rings 
with an expected line amplitude above a certain threshold are 
now discarded. 



Cross-correlation factor between the 4K line amplitudes 


Figure 7. Normalized histograms of the correlation coefficients 
of the 10 and 30 Hz 4-K line amplitudes. The amplitude is com¬ 
puted per ring and per bolometer from the two coefficients (sine 
and cosine) of a given line. For each pair of distinct bolometers 
(from 100 to 353 GHz), a correlation coefficient is computed be¬ 
tween the two amplitudes during the mission. The black normal¬ 
ized histogram shows the 10 Hz line correlation coefficients of 
the 903 (= 43 x 42/2) pairs. The red curve shows the 30Hz line 
histogram. The 30 Hz line is clearly correlated between different 
bolometers. This is the only line that shows a significant corre¬ 
lation. 


In contrast to the other lines, the 30 Hz line signal is corre¬ 
lated across bolometers (see Fig. 7). It is therefore likely that 
the 4-K line removal procedure leaves correlated residuals on 
the 30 Hz line. The consequence of this correlation is that the 
cross-power spectra between different detectors can show ex¬ 
cess noise at multipoles around £ ^ 1800^ (see the discussion in 
Section 1 of Planck Collaboration XVI 2014 and in Section 7 of 


The spacecraft spin period of one minute implies a correspondence 
between a TOI frequency of 30 Hz and C 1800. 
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Multipole I 

Figures. Temperature cross power-spectrum of the 217GHz 
detector sets 1 and 2 for the full mission (black) and yearly cuts 
(Year 1 in red and Year 2 in green), comparing the 2013 (dashed 
lines) and 2015 (continuous lines) data release. 


Planck Collaboration XV 2014). An example, computed with the 
anafast code in the HEALPix package (Gorski et al. 2005), is 
shown in Fig. 8. To mitigate this effect, we discard all 30 Hz res¬ 
onant rings for the 16 bolometers between 100 and 353 GHz for 
which the median average of the 30 Hz line amplitude is above 
10 aW. Doing so, the I - 1800 feature disappears. This issue 
is also addressed in Section 3.1 of Planck Collaboration XIII 
(2015). 

Figure 9 shows the fraction of discarded samples for each 
detector over the full mission. It gathers the flags at the sample 
level, which are mainly due to glitches and the depointing be¬ 
tween rings. It also shows the flags at the ring level, which are 
mostly due to the 4-K lines, but are also due to solar flares, big 
manoeuvres, and end-of-life calibration sequences, which are 
common to all detectors. The main difference from the nominal 
mission, presented in the 2013 papers, appears in the fifth survey, 
which is somewhat disjointed, due to solar flares arising with 
the increased solar activity, and to special calibration sequences. 
The full cold Planck-WFl mission lasted 885 days, excluding 
the calibration and performance verihcation (CPV) period of 
1.5 months. During this time, HFI data losses amount to 31 %, 
the majority of which comes from glitch flagging, as shown in 
Table 1. The fraction of samples flagged due to solar system ob¬ 
jects (SSO), jumps, and saturation (Planck Collaboration VIII 
2014) is below 0.1 %, and hence negligible. 

3 . 4 . Detector time response 

As noted in Planck Collaboration VII (2014) and Planck HFI 
Core Team (2011a), the detector time response is a key calibra¬ 
tion parameter for HFI. It describes the relation between the opti¬ 
cal signal incident on the detectors and the output of the readout 
electronics. This relation is characterized by a gain, and a time 
shift, dependent on the temporal frequency of the incoming op¬ 
tical signal. As in previous releases, it is described by a linear 
complex transfer function in the frequency domain, which we 


call the time transfer function. This time transfer function must 
be used to deconvolve the data in order to correct the frequency- 
dependent time shift, which otherwise significantly distorts the 
sky signal. The deconvolution also restores the frequency de¬ 
pendent gain. It is worth noting that: (i) the deconvolution sig¬ 
nificantly reduces the long tail of the scanning beam; (ii) it also 
symmetrizes the time response, which allows us to combine sur¬ 
veys obtained by scanning in opposite directions; and (iii) given 
that the gain decreases with frequency, the deconvolution boosts 
the noise at high frequency, as can be seen in Figs. 23 and 26. In 
order to avoid unacceptably high noise in the highest temporal 
frequencies, a phaseless low-pass filter is applied, with the same 
recipe as in Section 2.5 of Planck Collaboration VII (2014). This 
process results in a slightly rising noise in the high frequency 
part of the noise power spectrum, in particular for the slowest 
100 GHz detectors. This noise property is ignored in the map 
making process, which assumes white noise and low frequency 
noise. We note, however, that the 100 GHz bolometers data are 
not used for CMB analysis at smaller angular scales, due to their 
wider main beams. 

For this release, the time transfer function is based on the 
same model as the previous release. As a function of the angular 
frequency o), it is defined by 


TF(w) = S pha.se, Tstray), (3) 

where F{a>) is the term associated with the bolometer response, 
and //'(w; 3phase, Tstray) IS the analytic model of the electronics 
transfer function, whose detailed equations and parameters are 
given in Appendix A of Planck Collaboration VII (2014). The 
electronics term depends only on two parameters: the phase shift 
between the AC bias current and the readout sampling, encoded 
by the parameter 3phase; and the time constant Tso-ay associated 
with the stray capacitance of the cables connecting the bolome¬ 
ters to the bias capacitors (the stray capacitance enters in the 
h() term in the series of biters reported in table A. 1 of Planck 
Collaboration VII 2014). 

The bolometer time transfer function is described by the sum 
of bve single-pole low-pass functions, each with a time constant 
T, and an associated amplitude a,: 


F(co) = ^ 
(= 1,5 


aj 

1 -H icJTi 


(4) 


In this version, bve time constants are used rather than four, and 
the values of the parameters have been measured from a differ¬ 
ent combination of data than used previously. The extra low-pass 
function and the consequent parameter updates are motivated by 
the discovery of a time delay between the measured CMB dipole 
and the expected one. This time delay, or phase shift, is inter¬ 
preted as the effect of an extra low-pass function, not accounted 
for in the previous version of the deconvolution process. As we 
will see, some level of time delay in the dipole remains even after 
deconvolution of the bve low-pass functions. This small residual 
can generate problems in gain btting and dipole subtraction, and 
is treated more efficiently at the map making level (see Paper 2). 

In addition to the dipole shift, new information comes from 
the stacking of glitches induced by high energy particles hits. As 
described in Planck Collaboration X (2014), short glitches are 
due to direct interaction of particles with the bolometer grid or 
thermistor. The time response of short events is then representa¬ 
tive of the response to photons. Short glitches show a long tail, if 
enough events are stacked together. It was observed that the sig¬ 
nals obtained by stacking short glitches and by stacking Jupiter 
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Figure 9. Fraction of discarded data per bolometer (squares with black line). The fraction of data discarded from glitch flagging 
alone is shown with stars and the green line. The blue line with diamonds indicates the average fraction of discarded samples in 
valid rings. The two bolometers showing permanent random telegraphic signal (RTS), i.e. 143-8 and 545-3, are not shown, because 
they are not used in the data processing. 


Table 1. Overall budget of discarded data samples for the 50 valid bolometers. Two bolometers, one at 143 GHz and one at 545 GHz, 
cannot be used, due to a permanent random telegraphic signal (RTS) and are excluded from the statistics. A global average is given 
for glitches. Here, depointing denotes only the standard manoeuvres from one ring to another. Big manoeuvres are included in the 
common discarded rings (but not the 4-K line selection process). The RTS affects six bolometers episodically. The 4-K line selection 
process affects 20 bolometers (see Sect.3.3). The range of values obtained for different bolometers is given in the last column. Note 
that percentages do not add up to the total, since depointing, common rings, RTS and 4-K flagged samples are already flagged at 
the 20 % level due to glitches. 


Origin 

Mean fraction loss [%] 

Range [%] 

Glitch. 

. 20 

9-32 

Depointing . 

. 8 

8-8 

Common discarded rings . 

. 2 

2-2 

RTS . 

. 0 

0-4 

4-K . 

. 4 

0-16 

Total . 

. 31 

17^6 



Time [s] 


Figure 10. Impulse response of bolometer 143-6 to short glitches and to Jupiter. The plotted axis is linear within the range +10 ^ 
and logarithmic elsewhere. 


scans decay to zero with the same time constants (see Fig. 10). 
Since the short glitch stacking has the benefit of a much higher 
signal-to-noise ratio, it provides the most sensitive measurement 


of the longest time constants. Although the physical process of 
energy injection into the detector is different for microwave pho¬ 
tons and energetic particles, the heat dissipation is expected to 
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be the same. For this reason, it was decided to take only the time 
constants from the glitches, and not the associated amplitudes. 
The impact of the incomplete correction of the transfer function 
is discussed in Sect. 3.4.1. 

In summary, the values of the bolometer time transfer func¬ 
tion parameters, a, and t,, are measured with the following logic: 

- the two fastest time constants, ti and T 2 , and the associated 
amplitudes ai and 02 , are unchanged with respect to the pre¬ 
vious version, i.e., they are estimated from planet observa¬ 
tions; 

- the two longest time constants, T4 and T5, are estimated from 
short glitch stacking, together with the 05/04 ratio; 

- T3, 03, and 04 are fitted from Jupiter scans, keeping ti, oi, 
T2, and 02 fixed, while the value of 05 is set to keep the same 
ratio 05/04 as in the glitch data; 

- 05 is fitted from the CMB dipole time shift. It should be noted 
that the same dipole time shift can be obtained with different 
combinations of 05 and T5 and for this reason, T5 is recovered 
from the short glitches, and only 05 from the dipole time 
shift. 

This process was used for all channels from 100 to 217 GHz. 
For the 353 GHz detectors, different processing was needed to 
avoid strong non-optical asymmetries in the recovered scanning 
beam (see Appendix C for details). For the submillimetre chan¬ 
nels, 545 and 857 GHz, the time transfer function is identical to 
that of the 2013 release. The values of the time response param¬ 
eters are reported in Table C.l. The model and parameters will 
be improved by continuing this activity for future releases. 

3.4.1. Time response errors 

The beam model is built from time-ordered data deconvolved by 
the time transfer function. For this reason, and considering the 
constant rotation rate of Planck, the measured scanning beam 
absorbs, to a large extent, mismatches between the adopted time 
transfer function and the true one (Planck Collaboration VII 
2014). In this sense, errors and uncertainties in the time trans¬ 
fer function should not be propagated into an overall window 
function uncertainty, since the time response acts as an error- 
free regularization function. Biases and uncertainties are taken 
into account in the beam error budget. 

The beam window function correction of HFI’s angular 
power spectrum will not work on time scales longer than than 
roughly 0.5 s, given the size of the scanning beam map and 
the constant scan rate. Simulations with varying time response 
parameters show that errors at frequencies above 2 Hz are ab¬ 
sorbed in the scanning beam map and propagated correctly to 
the effective beam window function to better than the scanning 
beam errors. Errors on time scales longer than 0.5 s are not ab¬ 
sorbed by the scanning beam map, but propagate into the shifted 
dipole measurement and the relative calibration error budget (see 
Paper 2). 

Additionally, the time response error on time scales longer 
than 0.5 s can be checked by comparing the relative amplitude of 
the first acoustic peak of CMB anisotropies between frequency 
bands. The main calibration of HFI is performed with the CMB 
dipole (appearing in the TOI at 0.016 Hz), while the first acoustic 
peak at ^ ss 200 appears at 6 Hz. Table 3 in Planck Collaboration 
I (2015) shows that between 100 and 217 GHz the agreement is 
better than 0.3 %. 


3.4.2. Focal plane phase shift from fast Mars scan 

In December 2011, Planck underwent a series of HFI end-of-life 
tests. Among these, a speed-up test was performed increasing the 
spin rate to 1.4 rpm from the nominal value of 1 rpm. The test 
was executed on 7-16 December 2011, and included an obser¬ 
vation of Mars. Right after the test, a second Mars observation 
followed, at nominal speed. The main result of this test was the 
ability to set the real position of the detectors in the focal plane, 
which, when scanning at constant speed, is completely degener¬ 
ate with a time shift between bolometer data and pointing data. 
This time shift was supposed to be zero. During the test, it was 
found that the deconvolved bolometer signals peaked at a dif¬ 
ferent sky position for the nominal scans and for the fast scans. 
This discrepancy was solved by introducing a time-shift between 
bolometer readout and pointing data. This time-shift is not the 
same for all detectors, and ranges from 9.5 to 12.5 ms, with a 
corresponding position shift ranging from 3.'3 to 4.'5 in the scan¬ 
ning direction. Due to time constraints, it was not possible to 
observe Mars twice with the full focal plane, so it was decided 
to favour the CMB channels and the planet was observed with 
all the 100, 143, 217, and 353 GHz detectors. For the other de¬ 
tectors, we used the average time-shift of all the measured detec¬ 
tors. Notably, this detector-by-detector shift resulted in a better 
agreement of the position in the focal plane of the two PSBs in 
the same horn. 

Working with fast spin-rate deconvolved data is complicated 
by aliasing effects. For this reason the fitting procedure followed 
a forward sense approach, by modelling the signal with a beam 
centred in the nominal position, then convolving the fast and 
nominal timelines with the time transfer function, and fitting the 
correct beam centres by comparing data and model for both fast 
and nominal spin-rates. For the same reason, a direct compar¬ 
ison of deconvolved timeline has not proven to provide better 
constraints on the time-response parameters. 

4. Planets and main beam description 

We follow the nomenclature of Planck Collaboration VII (2014), 
where the “scanning beam” is defined as the coupled response of 
the optical system, the deconvolved time response function, and 
the software low-pass filter applied to the data. The “effective 
bea”m represents the averaging of signal due to the scanning 
of the telescope and mapmaking, and varies from pixel to pixel 
across the sky. 

Here we redefine the “main beam” to be the scanning beam 
out to 100' from the beam axis. The sidelobe structure at this 
radius is dominated by diffraction at the mirror edges and falls 
as oc 0^^, where 0 is the angle to the main beam axis. The main 
beam is used to compute the effective beam and the effective 
beam window function, which describes the filtering of sky sig¬ 
nals. 

The smearing of the main beam cannot be significantly re¬ 
duced without boosting the high-frequency noise. The regular¬ 
ization function (a low-pass filter) chosen has approximately the 
same width as the instrumental transfer function. The deconvo¬ 
lution significantly reduces the long tail of the scanning beam 
(see Fig. 11). The deconvolution also produces a more symmet¬ 
ric time response, so residual “streaking” appears both ahead and 
behind the main beam, though ahead of the beam it is at a level 
of less than 10 of the peak response. 

The “far side lobes” are defined as the response from 0 > 5°, 
roughly the minimum in the optical response. The response be- 
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Figure 11. A scan across Jupiter with bolometer 143-lb, pre-deconvolution in black, post-deconvolution in red. 


gins to rise as a function of angle 9 beyond this due to spillover. 
The far sidelobes are handled separately from the beam effects 
(see Sect. 4.6 for justification and Paper 2 for details). 

Observations of planets are used to estimate the main beams 
and to calibrate the 545 and the 857 GHz channels. The main 
beam is needed to correct for the filtering of the CMB sky 
by the instrument, details of which can be found in Planck 
Collaboration VII (2014). There were several key changes in the 
reconstruction of the main beam since the 2013 data release, de¬ 
scribed in detail in Appendix B. 

- The TOI from Saturn and Jupiter observations are merged 
prior to B-spline decomposition, taking into account resid¬ 
ual pointing errors and variable seasonal brightness. This is 
achieved by determining a scaling factor and a pointing off¬ 
set by fitting the TOI to a template from a previous estimate 
of the scanning beams. We iterate the planet data treatment, 
updating the template with the reconstructed scanning beam. 
The process converges in five iterations to an accurcay of 
better than 0.1 % in the effective beam window function. 

- Steep gradients in the signal close to the planet reduce the 
completeness of the standard glitch detection and subtraction 
procedure, so the planet timelines are deglitched a second 
time. 

- The beam pipeline destripes the planet data, estimating a sin¬ 
gle baseline between 3° and 5° before the peak for each scan¬ 
ning circle. Baseline values are smoothed with a sliding win¬ 
dow of 40 circles. The entire scanning circle is removed from 
the beam reconstruction if the statistic in the timeline region 
used to estimate the baseline is far from Gaussian. 

- The main beam is now recontructed on a square grid that 
extends to a radius of 100' from the centroid, as opposed to 
40' in the 2013 data release. The cutoff of 100' was chosen 
so that a diffraction model of the beam at large angles from 
the centroid predicts that less than 5 x 10“^ of the total solid 
angle is missing. 

- The scanning beam is constructed by combining data from 
Saturn observations, Jupiter observations, and physical op¬ 
tics models using GRASP software. 

- No apodization is applied to the scanning beam map. 

The update of the time response deconvolved data has 
slightly changed the scanning beam, the effective beam solid an¬ 
gles, and the effective beam window functions. 


Table 2. Band-average scanning beam solid angle (Gsb) and 
Monte Carlo-derived errors (AGmc) including noise, residual 
glitches, and pointing uncertainty. 


Band fJsB 

[GHz] [arcmin^] 


100 . 104.62 0.13% 

143 . 58.80 0.07% 

217 . 26.92 0.13% 

353 . 25.93 0.09 % 

545 . 25.23 0.08 % 

857 . 23.04 0.08 % 


4. 1. Hybrid beam model 

A portion of the near sidelobes was not accounted for in 
the effective beam window function of the 2013 data (Planck 
Collaboration VII 2014; Planck Collaboration XXXI 2014). To 
remedy this, the domain of the main beam reconstruction has 
been extended to 100', with no apodization. Saturn data are used 
where they are signal-dominated. Where the signal-to-noise ratio 
of the Saturn data falls below 9, azimuthally binned Jupiter data 
are used. At larger angles, below the noise floor of the Jupiter 
data, we use a power law (oc 6^^), whose exponent is derived 
from GRASP, to extend the beam model to 100'. Figure 12 shows 
a diagram of the regions handled differently in the hybrid beam 
model. A summary of the solid angles of the hybrid beams is 
shown in Table 2. Figure 13 shows a contour plot of all the scan¬ 
ning beams referenced to the centre of the focal plane. 

4.2. Effective beams and window functions 

As described in Mitra et al. (2011) and Planck Collaboration VII 
(2014), the FEBeCoP code is used to compute the effective beam 
(the scanning beam averaged over the scanning history) and the 
FEBeCoP and Quickbeam codes are used to compute the effec¬ 
tive beam window functions. Statistics of the effective beams are 
shown in Table 3. 

Using Saturn and Jupiter to reconstruct the main beam in¬ 
troduces a small bias due to the large disc size of the plan¬ 
ets as compared to Mars (see Fig. 8 of Planck Collaboration 
VII 2014). Additionally, Saturn’s ring system introduces a slight 
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Figure 12. Scanning beam map for detector 143-6 with a rough 
illustration of the regions that are handled with a different data 
selection or binning. 


4.3. Beam error budget 

As in the 2013 release, the beam error budget is based on an 
eigenmode decomposition of the scatter in simulated planet ob¬ 
servations. A reconstruction bias is estimated from the ensemble 
average of the simulations. We generate 100 simulations for each 
planet observation that include pointing uncertainty, cosmic ray 
glitches, and the measured noise spectrum. Simulated glitches 
are injected into the timeline with the correct energy spectrum 
and rate (Planck Collaboration X 2014) and are detected and 
removed using the deglitch algorithm. Noise realizations are de¬ 
rived as shown in Sect. 5.3.1. Pointing uncertainty is simulated 
by randomizing the pointing by 1!'5 rms in each direction. 

The improved signal-to-noise ratio compared to 2013 leads 
to smaller error bars; for instance, at £ - 1000 the uncertainties 
on are now (2.2,0.84,0.81) x lO^^ for 100, 143, and 217 GHz 
frequency-averaged maps respectively, reduced from the previ¬ 
ous uncertainties of (61,23,20) x lO^"^. 

A “Reduced Instrument Model” (RIMO; see Appendix A) 
containing the effective B{{) for temperature and polarization de¬ 
tector sets, for auto- and cross-spectra at 100 to 217 GHz, is in¬ 
cluded in the release, for a sky fraction of 100 %. Another RIMO 
is provided for a sky fraction of 75 %. They both contain the first 
five beam error eigenmodes and their covariance matrix, for the 
multipole ranges [0, ^max] with = 2000, 3000, 3000 at 100, 
143, and 217 GHz respectively (instead of 2500, 3000,4000 pre¬ 
viously). These new ranges bracket more closely the ones ex¬ 
pected to be used in the likelihood analyses, and ensure a better 
determination of the leading modes on the customized ranges. 

As described in Appendix A7 of Planck Collaboration XV 
(2014), these beam window function uncertainty eigenmodes are 
used to build the C({) covariance matrix used in the high-^ an¬ 
gular power spectrum likelihood analysis. It was found that the 
beam errors are negligible compared to the other sources of un¬ 
certainty and have no noticeable impact on the values or associ¬ 
ated errors of the cosmological parameters. 


frequency-band dependence of the effective size of the planetary 
disc (Planck Collaboration XXXIV 2015). Because of Saturn’s 
small size relative to the HFI beams, for { < 4000 the symmet¬ 
ric part of Saturn’s shape dominates the window function. There 
is a small variation with HFI band in the apparent mean size 
of Saturn due to the different ring system temperatures, ranging 
from 9725 at 100 GHz to 1072 at 857 GHz. Because this varia¬ 
tion introduces a bias of less than 2x10“^ in for £ < 4000, we 
ignore it and use a correction derived for a mean 9!'5 Saturn disc 
for all bands. 

The two effective-beam codes handle temperature-to- 
polarization leakage in slightly different ways. The dominant 
leakage comes from differences in the scanning beams of the 
polarization-sensitive detectors (See Sect. 4.7). The Quickbeam 
code mimics the data model of the mapmaking code and assumes 
that every polarization-sensitive detector at a given frequency 
has the same beam shape, and thus produces a single effective 
beam window function. Temperature-to-polarization leakage is 
handled later as a set of parameters in the likelihood of the angu¬ 
lar power spectra. The FEBeCoP code produces effective beam 
window functions for the polarized power spectra that account 
for differences in the main beam. However, these are computed 
as the average power leakage of a given sky signal from temper¬ 
ature to polarization. Hence, these polarized window functions 
are not strictly instrumental parameters, since they rely on an 
assumed fiducial temperature angular power spectrum. 


4.4. Consistency of beam reconstruction 

To evaluate the accuracy and consistency of the beam recon¬ 
struction method, we have compared beams reconstructed using 
Mars for the main beam part instead of Saturn, and using data 
from Year 1 or Year 2 only. We compared the window functions 
obtained with these new beams to the reference ones; new Monte 
Carlo simulations were created to evaluate the corresponding er¬ 
ror bars, and the results are shown in Figs. 14 and 15. Note that 
the effective beam window functions shown in these hgures are 
not exactly weighted for the scan strategy, rather we plot them in 
the raster scan limit, where B^ is defined as the sum over m the 
B^^. The eigenmodes were evaluated for these different data sets 
allowing us to compare them with the reference beam using ax^ 
analysis. The discrepancy between each data set and the refer¬ 
ence beam is fitted using Vdof=5 eigenmodes using fmax =1200, 
2000, and 2500, for the 100, 143, and 217 GHz channels, respec¬ 
tively. Thex^ is then defined as: 

A^dof 

= ^(c,/4,)2 (5) 

1=1 

where c, is the fit coefficient for eigenvector i, with eigenvalue 
/!, . For each data set, the p-value to exceed this;^f^ value for ax^ 
distribution with Vdof degrees of freedom is indicated in Table 4. 
A high p-value indicates that the given data set is consistent 
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Cross-scan [°] 


Figure 13. B-spline hybrid scanning beams reconstructed from Mars, Saturn, and Jupiter. The beams are plotted in logarithmic 
contours of -3, -10, -20, and -30 dB from the peak. PSB pairs are indicated with the a bolometer in black and the b bolometer in 
blue. 


with the reference beam within the simulation-determined error 
eigenvectors, with 100 % indicating perfect agreement. 

We find excellent agreement between the yearly and nominal 
beams for the SWB bolometer channels. However, in order to 
obtain reasonable beam agreement for the PSB detector sets, we 
find that the yearly beam errors must be scaled by a factor of 
4. We conclude that there is an unknown systematic error in the 
beam reconstruction that is not accounted for in the simulations 
used to estimate the Monte Carlo error bars. We assign a scaling 
factor of 4 to the error eigenmodes to account for this systematic 


uncertainty. The additional error scaling has a negligible effect 
on the cosmological parameters. 

4.5. Colour correction of the beam shape 

In general, the measured beam shape is a function of the spectral 
energy distribution (SED) of the measurement source. This is of 
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Table 3. Mean values of effective beam parameters for each HFI frequency. The error in the solid angle ctq comes from the scanning 
beam error budget. The spatial variation AQ is the rms variation of the solid angle across the sky. The reported FWHM is that of the 
Gaussian whose solid angle is equivalent to that of the mean effective beam. Qi and Q 2 are the solid angles contained within circles 
of radius 1 and 2 FWHM, respectively (used for aperture photometry as described in Appendix A of Planck Collaboration XXVIII 
2014). 


Band ctq AQ FWHM Q 2 

[GHz] [arcmin^] [arcmin^] [arcmin^] [arcmin] [arcmin^] [arcmin^] 


100 . 106.22 0.14 0.20 9.69 100.78 106.03 

143 . 60.44 0.04 0.20 7.30 56.97 60.21 

217 . 28.57 0.04 0.19 5.02 26.46 28.46 

353 . 27.69 0.02 0.20 4.94 25.32 27.53 

545 . 26.44 0.02 0.21 4.83 24.06 26.09 

857 . 24.37 0.02 0.12 4.64 22.58 23.93 



MurnpoO' Mukpol*' MuDipol* > MuHIpoK > UuhfiOlt < Multool*' 

Figure 14. Comparison of Year 1 and Year 2 based beams with the reference beam (Full Mission). Window functions are calculated 
using Quickbeam in raster scan conhguration. Error bars are computed using MC simulations for Year 1 and Year 2, taking the 
maximal one for each. 




1(7 

3(7 


Multipole ^ 


Multipole ^ 



Figure 15. Comparison of Mars-based beams with the reference 
beam (Full Mission). Mars-based beams are constructed from 
Mars data for the main beam (< 10') and Saturn and Jupiter 
data for the larger scales. Window functions are calculated using 
Quickbeam in raster scan conhguration. Error bars are computed 
using MC simulations for Mars. 


Table 4. P-values in percent for the ;^^^comparison of the nomi¬ 
nal beam with Year 1 and Year 2 data sets, following the dehni- 
tion in Eq. (5). 


Band 

Year 1 

Year 2 

100-dsl _ 

. 66.95 

92.44 

100-ds2 _ 

. 85.12 

75.99 

143-dsl _ 

. 72.24 

99.83 

143-ds2 _ 

1.14 

16.81 

143-5. 

. 94.88 

97.57 

143-6. 

. 96.14 

98.78 

143-7. 

. 93.78 

95.67 

217-dsl _ 

. 78.12 

69.74 

217-ds2 _ 

. 27.33 

30.87 

217-4 . 

. 95.99 

68.44 

217-1 . 

. 99.08 

97.49 

217-2 . 

. 97.18 

98.86 

217-3. 

. 97.59 

99.07 


particular concern, because we measure the beam on a source 

with a roughly Rayleigh-Jeans SED, yet we use the effective Planck Collaboration VII (2014) described possible levels of 
beam window function to correct the CMB. this bias derived using GRASP calculations with the pre-launch 
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telescope model (Maffei et al. 2010; Tauber et al. 2010). The 
pre-launch calculations did not agree well enough with the data 
to allow a direct application of the colour correction. A new 
telescope model based on flight data, and presented in a forth¬ 
coming paper, predicts beams that agree better with the data 
at 100-217 GHz, but show worse agreement at 353 GHz. The 
new model predicts a different beam shape colour correction, 
though at a similar order of magnitude to that shown in Planck 
Collaboration VII (2014). This work is ongoing and will be com¬ 
pleted after the 2015 release. 

We search for a possible signature of the colour-correction 
effect in the data. First, we check for consistency in the CMB 
angular power spectra derived from different detectors and fre¬ 
quency bands using the SMICA algorithm, in order to find dis¬ 
crepancies in beam shape and relative calibration of different 
data subsets using the CMB anisotropies (Planck Collaboration 
XI2015). There are hints of differences that are orthogonal to the 
beam error eigenmodes and appear as changes in relative cali¬ 
bration, but the preferred changes are at the 0.1% level. Because 
this correction is of the same order of magnitude as the relative 
calibration uncertainty between detectors, we treat it as insignif¬ 
icant. 

Second, we look at the relative calibration between detectors 
within a band for compact sources as a function of the source 
SED. In this method, any discrepancies due to solid angle vari¬ 
ation with SED are degenerate with an error in the underlying 
bandpass. Additionally, because of the lack of bright, compact 
sources with red spectra, this method is limited in its signal-to- 
noise ratio. In the limit that any detected discrepancy is com¬ 
pletely due to solid-angle variation with colour, at the level of 
1 % in solid angle we do not detect variations consistent with the 
physical optics predictions. 

Given the lack of measurement of this effect at the current 
levels of uncertainty, as well as uncertainties in the modelling of 
the telescope, we note that this effect may be present at a small 
level in the data, but we do not attempt to correct for it or include 
it in the error budget. 

4.6. Effective beam window function at iarge anguiar scaies 

The main beam derived from planet observations and GRASP 
modelling extends 100' from the beam axis, so the effective 
beam window function does not correct the filtering of the sky 
signal on larger scales (approximately multipoles { < 50). In 
practice, due to reflector and baffle spillover, the optical response 
of HEI extends across the entire sky (Tauber et al. 2010; Planck 
Collaboration XIV 2014). According to GRASP calculations, the 
far-sidelobe pattern (ESL) further than 5° from the beam cen¬ 
troid constitutes between 0.05 % and 0.3 % of the total solid an¬ 
gle. To first order, this is entirely described by a correction to 
the calibration for angular scales smaller than the dipole (Planck 
Collaboration XXXI 2014). 

More correctly, the far-sidelobe beam (defined here as the 
optical response more than 5° from the beam axis) filters large- 
angular-scale sky signals in a way that is coupled with the scan 
history of the spacecraft, and given a perfect measurement of 
the far-sidelobe beam shape, an effective beam window func¬ 
tion could be constructed for the large angular scale with the 
same procedures used for the main beam’s effective beam win¬ 
dow function. 

We rely entirely on GRASP calculations to determine the 
large-scale response of the instrument. These calculations are fit¬ 
ted to survey difference maps (see Planck Collaboration XXXIV 
(2015)), where the single free parameter is the amplitude of the 


GRASP far-sidelobe model beam. Because of a combination of 
the low signal-to-noise ratio of the data and uncertainty in the 
GRASP model, the errors on the fitted amplitude are of the or¬ 
der of 100 %. While ground-based measurements have given us 
confidence in our limits on total spillover, we note that there 
is a large uncertainty in the exact shape of the far sidelobes. 
Concerning a possible evolution of the far sidelobe during the 
mission, we note that the background, as measured by the total 
power on the bolometers, decreases then stabilizes as a function 
of time. This proves that dust deposition on the Planck mirrors is 
negligible. Hence, the far sidelobes are likely to be very stable. 

Pig. 16 shows the correction to the effective beam window 
function due to the best-fit far-sidelobe model. The error bars are 
from the amplitude error in the fit (Planck Collaboration XXXIV 
2015) and they show that the uncertainty is larger than the pre¬ 
dicted beam window function correction except at the dipole 
scales {i - 1), aside from a marginally significant, but negligibly 
small, bump at f = 5. The filtering of the CMB signal with the 
far-sidelobe beam is thus insignificant and can be ignored as a 
component of the effective beam window function. A map-level 
correction of Galactic contamination pickup in the far-sidelobes 
is more appropriate. However, given the large uncertainties in the 
far-sidelobe model shape, we choose to neglect this correction, 
and note that residuals from Galactic pickup may be present in 
the map. The only far sidelobe effect that is corrected in the data 
is an overall calibration correction, as described in Paper 2. 



Figure 16. Estimate of the effective beam window function cor¬ 
rections due to far-sidelobe response (more than 5° from the 
beam axis). Shading indicates +lcr errors. The curves in the 
353 GHz panel are divided by a factor of 10; for this frequency 
the far sidelobe response amplitude is very low compared to the 
estimated errors. 


4.7. Cross-poiar response and temperature-to-poiarization 
ieakage 

The data model used for the HEI polarization reconstruction 
assumes that the entire cross-polar response of each PSB is 
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Figure 17. A simulation of the ratio of temperature-to-polarization leakage, given the main beam mismatch, compared to EE and 
BB signal. 


due to the detector itself, and so the beam shape of the cross- 
polar response is exactly the same as the co-polar beam shape. 
The corrugated feed horns exhibit some internal cross-polar 
response (Maffei et al. 2010). GRASP simulations show that 
this cross-polar response is at the level of 0.1-0.5%, consid¬ 
erably smaller than the detector cross-polar response (2-5%). 
Simulations using FEBeCoP show that the effect of ignoring 
the cross-polar optical beam shape results in a smoothing of 
polarization power spectra equivalent to an additional 10-20" 
Gaussian, which we neglect. Additionally the data model as¬ 
sumes identical beams for every PSB used to reconstruct the 
polarization; this creates some temperature-to-polarization leak¬ 
age. Main-beam leakage dominates over the differential time re¬ 
sponse tails. Figure 17 shows a simulation of CMB temperature 
leakage into polarization power spectra using FEBeCoP, given 
the measured scanning beams within a band. 

5. Validation and consistency tests 

Here we describe some of the tests that have been done to vali¬ 
date the quality of the cleaned TOI. Other tests at the map level 
are described in Paper 2. First we discuss the impact of the ADC 
correction. Then we analyse how each step of the TOI process¬ 
ing pipeline can alter the resulting power spectra and also study 
the noise properties. Finally we estimate the filtering function of 
the TOI pipeline with end-to-end simulations. 

5.1. ADC residuals 

5.1.1. Gain consistency at ring level 

Using the undeconvolved TOI at the ring level, we can monitor 
the quality of the ADC correction with respect to the stability of 
the gain. For that purpose we can measure the relative gain of 
parity plus (g+) and parity minus (g_) samples (alternating sam¬ 


ples). The two parities sample a very different part of the ADC 
scale, so that a gain mismatch is a diagnostic of the ADC nonlin¬ 
earity correction. During a ring, the sky signal from both parities 
should be almost identical. We thus correlate a phase-binned ring 
(PBR) made of parity plus with the average PBR (made of both 
parities) and obtain g+. We similarly obtain g_. The gain half¬ 
difference (g+ - g-)/2 is shown as a function of the ring num¬ 
ber in Fig. 18 for a representative selection of four bolometers. 
The improvement obtained with ADC correction is significant 
with a root-mean-square dispersion decreasing by a factor of 2 
to 3. Only a handful of bolometers show some discrepancies at 
the 10“^ level after the ADC correction, namely 143-3b, 217-5b, 
217-7a, 217-8a, and 353-3a. This internal consistency test at the 
ring level is not sensitive to gain errors common to both parities. 
However, it agrees qualitatively with the overall Bogopix gain 
variations shown in Paper 2. 

5.1.2. Parity map spectra 

Two independent sky maps are produced from the two TOI pari¬ 
ties; the differences between these are again expected to capture 
the ADC residual effects. Figure 19 compares full-sky spectra 
derived from maps built on raw and ADC-corrected data. Glitch- 
flagged samples have been removed. The behaviour in Xjt} at 
low I is similar to the behaviour of simulated data as shown in 
Sect. 2.5. 

These low-^ residuals, although much improved with respect 
to the 2013 delivery, are not yet fully under control. At the time 
of the 2015 delivery, work is continung to estimate how this ef¬ 
fect propagates into the cross-polarization spectra, and to give 
a coherent description of the different low-^ systematics still 
present in HFI 2015 data. For this reason, among others, it has 
been decided to postpone the release of HFI low t polarized data. 
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Figure 18. The relative gain difference between parity + and - samples is shown as a function of the ring number. A boxcar average 
of 101 rings has been applied. The purple data points were processed with the 2013 pipeline, without ADC correction, while the red 
points refer to the present data release with ADC correction. The median error bars are of the order of 10“^ for the ADC-corrected 


TOI. 


5.2. Filtering effects 

It is customary to evaluate the effects produced by the pipeline 
on the signal as a filtering function, depending on the considered 
angular scale. Here we evaluate whether or how each pipeline 
module (see Fig. 1) filters the signal. 

- ADC correction. The ADC nonlinearities add a spurious 
component to the power spectra. Correcting for it removes 
this component, reducing the C/ by an amount smaller than 
1 %, that depends on the channel, and is very flat in £. 

- Glitch removal/flagging. Using simulations, Planck 
Collaboration X (2014) have shown that there is no sig¬ 
nificant bias (< 10 induced on the signal by the glitch 
removal procedure, thanks to the joint estimation of glitch 
tail and sky signal. The sky signal is estimated at the ring 
level using spline interpolation in order to correct for signal 
variations within ring pixels. However, the pointing jitter, 
combined with a strong signal gradient, can influence the 
glitch detection rate, mostly around the Galactic plane and 
planets. We have made simulations in order to evaluate an 
equivalent filtering effect using the glitch module called 
despike (Planck Collaboration X 2014). For CMB chan¬ 


nels, the impact is clearly negligible. For submillimetre 
channels, the signal (mostly the interstellar dust emission) is 
changed at the level of 10 

- AC modulation baseline subtraction. This has no impact on 
the signal, since it affects the TOI at the modulation fre¬ 
quency, which is cut out at the end of the TOI processing. 

- Thermal decorrelation. This involves subtracting a TOI that 
is filtered at the minute timescale. The computation of the 
offsets per ring (in the destriping during the map making pro¬ 
cess) absorbs the long-term variability (typically longer than 
an hour) but not the mid-term variability (between 1 min and 
1 h). The impact of the thermal decorrelation is to improve 
the accuracy of the destriping offsets. It is therefore not a fil¬ 
tering effect. The maximal impact of the thermal decorrela¬ 
tion is estimated by comparing Cf using data processed with 
and without the thermal decorrelation. In Fig. 20, we see that 
the cross-spectrum of data that are fully decorrelated in time 
(as in Year 1 x Year 2) shows differences less than 10 % for 
£ < 2000. The very-high-f’ behaviour reflects the noise level 
of the reference spectrum. Consistency tests in the likelihood 
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Multi pole I 


Figure 19. Half-parity map difference spectra for the same bolometers as in Fig. 18. Blue: 2013 release without ADC correction. 
Red: the present data release with ADC correction. Below ( - 100 the noise should be flat when there is no ADC bias. 



Figure 20. Relative differences of the C(, computed on 60% 
of the sky, using data processed with and without the thermal 
decorrelation. Note that the zero-level has been recomputed in 
all cases to allow the comparison to be made. The amplitudes of 
the negative values are shown by dotted lines. 


paper (Planck Collaboration XI 2015) are powerful enough 
to assess that the thermal residuals are negligible. 


- 4-K line removal. The maximal impact is estimated by com¬ 
paring Ce using data analysed with and without the removal 
of the 4-K line frequencies (keeping the same dehnition of 
valid rings in both cases). In Fig. 21 we see that only two 
lines have a strong impact: the 10 Hz line influences the 
power spectrum at f 600 and the 30 Hz ai i ^ 1800. 
For the 10 Hz line, the impact on the power spectrum is only 
seen in the auto-spectrum while the 30 Hz line remains in 
cross-spectra between two bolometers or between half rings. 
When the cross-spectrum is computed using data that are 
fully decorrelated in time (as in Year 1 x Year 2), no residual 
is seen. Consistency tests in the likelihood paper are pow¬ 
erful for assessing that the 4-K line residuals are negligi¬ 
ble, thanks to the extra ring discarding scheme described in 
Sect. 3.3. 

- Nonlinearity correction. The impact of this correction on 
measurements of the CMB has been propagated through the 
whole pipeline. Figure 22 shows that the angular power spec¬ 
trum is underestimated by less than 0.3 % between ( - 20 
and { - 2000 if no correction is applied, while still keeping 
the same calibration coefficients. We estimate that the non¬ 
linearity correction is accurate within a 10% uncertatinty. 
Thus, the uncertainty of the bolometer nonlinearity correc- 
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Figure 21. Relative differences of the Q, computed on 60 % of 
the sky, using data analysed with and without removal of the 4-K 
line frequencies. 


tion translates into a hltering function uncertainty below the 
10 “^ level over all angular scales. 



Figure 22. Relative differences of the Q, computed on 60 % of 
the sky, using data processed with and without the nonlinearity 
correction of the bolometer response. The example shown here 
is for bolometer 217-1. 


- The Fourier transform module does of course filter the TOI. 
It deconvolves the time transfer function (Sect. 3.4) but also 
low-pass biters the data (Planck Collaboration VII 2014). It 
is important for high temporal frequencies, but affects all fre¬ 
quencies at the 10 level down to the spin frequency where 
the biter is set to 1. The bitering effect is captured in the ef¬ 
fective beam for scales of up to 1-2°. However, there could 
remain an effective bitering effect for larger angular scales. 

In order to test the global pipeline bitering effect that could be 
induced by the interaction of several systematic effects, we resort 
to complete simulations, described in detail in Sect. 5.4. 


5.3. Noise analyses 

5.3.1. Noise Spectrum 

Noise spectra are best studied at the ring level using the re¬ 
dundancy built into the scanning strategy. The method is de¬ 
scribed in Planck Collaboration VI (2014) and consists of re¬ 
moving the sky-signal ring average from the TOI for each ring. 
Figure 23 shows the amplitude spectral density for representa¬ 
tive HFI bolometers. This is modelled in Sect. 5.3.3 for the pur¬ 
pose of providing accurate inputs for simulations. For bolome¬ 
ters with relatively long time constants (especially at 100 GHz), 
the electronic noise is boosted at high frequency. 

5.3.2. Noise stationarity 

The noise stationarity is best studied by measuring the rms de¬ 
viation of a TOI sample at the ring level, at the end of the TOI 
processing, where the average ring signal has been subtracted. 
The statistics are computed for samples which are valid, far from 
the Galactic plane, and far from bright point sources, in order 
to avoid strong residual gradients (see Planck Collaboration ES 
(2013)). The bias induced by the signal subtraction is corrected 
for; the average value (after calibration) over the mission is given 
in Table 5. This is compatible with the values quoted in Table 4 
of Planck Collaboration VI (2014). 

Table 5. HFI average total noise. The rms noise is given for one 
bolometer. The averaging between bolometers of the same band 
is done with the same weighting scheme as used in the mapmak¬ 
ing procedure. The band average (4th column) refers to the sen¬ 
sitivity of the collection of Nh bolometers at the same frequency. 


Band 

Vb 

Bol. total noise 

Band total noise 

Units 

100 

8 

113.3 

40.0 

^Kcmb 

143 

11 

57.5 

17.3 

pKcMB 

217 

12 

83.2 

24.0 

pKcMB 

353 

12 

282.0 

81.4 

^Kcmb 

545 

3 

45.5 

26.3 

kJy sr-‘ s‘/2 

857 

4 

49.2 

24.6 

kjy sr“* s'/^ 


Figure 24 shows trends of the total noise versus ring num¬ 
ber before (2013 release) and after the ADC correction (this re¬ 
lease), and histograms of the noise. In most cases, there is a sig- 
nihcant decrease in the relative width of the histogram when the 
ADC correction is included, though that decrease is more signif¬ 
icant for the low frequency channels, and since the distributions 
are non-Gaussian, very large decreases (or increases) indicate a 
large change in the structure of the distribution. Nevertheless, 
there are several cases where a jump in the noise is removed by 
the ADC correction (143-2a, 143-4a, 353-2). Other noise jumps 
(mostly at the sub-percent level except for bolometer 353-3a) re¬ 
main unexplained. The noise of some bolometers shows a linear 
trend with time, but at the sub-percent level. For most bolome¬ 
ters, the noise is stationary below the percent level. 

We have investigated the origin of the jumps in the noise 
level (e.g., for bolometer 353-2 at ring 5000 in the lower panel 
of Fig. 24). Most of the cases are not due to the ADC correc- 
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Figure 23. Calibrated noise spectrum for one bolometer in each HFI frequency channel. The spectrum is the average of all valid 
ring spectra of the second survey. The comb-like structure of absorption features occurs at multiples of the spacecraft spin frequency 
(16.7mHz), and is due to the signal removal procedure. At lOHz and above there are some residuals of the 4-K cooler lines only 
associated with this particular computation. 


tions (an exception is the 143-4a bolometer in Fig. 24). The extra 
noise component has been analysed by computing power spec¬ 
tra for selected rings. This falls into two categories: (1) for some 
detectors including 353-3a (the strongest case), it corresponds 
to an increase of the white noise component; (2) for other de¬ 
tectors, the extra component in the power spectrum is a bump 
concentrated in the 0.1 to 1 Hz range, which could be a very low 
level of RTS not detected by the non-Gaussianity tests that are 
routinely performed (Planck Collaboration VI 2014). Note that 
the small changes of noise levels with time are not taken into ac¬ 
count by the mapmaking process, because the loss of optimality 
can be neglected. 

5.3.3. Noise modelling 

Individual ring noise Power Specral Densities (PSDs) suffer 
from realization noise and signal residuals. For the purpose of 
simulating realistic noise from the PSDs, they need to be regu¬ 
larized by fitting a smooth model to the noisy PSD. To this end, 
we have constructed a physically-motivated model of the noise 
PSD consisting of 

1 . a photon/phonon noise component that is subject to the bolo- 
metric transfer function (suppressed above 10 Hz); 

2 . an electronic noise component that is only subject to the digi¬ 
tal low-pass filtering (abrupt suppression close to the Nyquist 
frequency); 

3. a 1// slope comprising the low-frequency thermal fluctua¬ 
tions, glitch residuals, and other low frequency fluctuations. 


The resulting noise model has four degrees of freedom: three 
amplitudes and one slope. 

In addition to the noise autospectrum components, we de¬ 
tect the presence of a correlated noise component, i.e., a com¬ 
mon mode between two polarization-sensitive bolometers in the 
same horn (e.g., 100-la and 100-lb). The common mode was 
measured by analysing the noise of sum, {a + b), and difference, 
(a - b), time streams. If the bolometers were entirely indepen¬ 
dent, the two composite streams would have equivalent noise 
power. Instead, we systematically find that (a + b) has a higher 
noise than (a - b). The measurements for 100 GHz are presented 
in Fig. 25. There are indications that similar common modes ex¬ 
ist across the focal plane, but for the purpose of noise simulation, 
the ajb correlated noise that directly affects polarization sensi¬ 
tivity is the most important one. Examples of noise model fits 
are presented for two bolometers in Fig. 26. Note that a common 
mode tracing the global lOOmK temperature fluctuations at the 
1 min time scale has already been removed from the processed 
TOI, by reference to the dark bolometers. 

5.4. End-to-end simulations 

End-to-end simulations are created by feeding the TOI process¬ 
ing pipeline with simulated TOI to evaluate and characterize 
its overall transfer function and the respective contribution of 
each individual effect. Simulated TOIs are produced by apply¬ 
ing the real mission scanning strategy to a realistic input sky 
specified by the Planck Sky Model (PSM vl.7.7, Delabrouille 
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Figure 24. Noise stationarity for a selection of four bolometers. The left panels show the total noise for each bolometer (dots) as 
a function of ring number. The total noise is the white noise spectral density giving the same rms per sample as the measured one. 
The solid line shows a running box average. The black dots are from the 2013 data release and the blue dots from this release. The 
right panels show histograms of the left trends. The box gives the width of the distribution at half maximum, as measured on the 
histogram, normalized to the mean noise level. The time response deconvolution has changed between the two data releases, so the 
absolute noise level is different. 
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Figure 25. Example of measured summed, (a + b), time stream 
noise power spectral densities (PSDs; dashed lines) and the dif¬ 
ference in the PSD between (a+b) and (a-b) (solid lines). These 
spectra are bin-by-bin median values across the entire mission 
and were derived by first measuring the sample autocovariance 
as described in Planck Collaboration XII (2015). 


et al. 2013) containing a lensed CMB realization, Galactic dif¬ 
fuse foregrounds, and the dipole components (but without point 
sources). To this scanning TOI, we then add a white-noise com¬ 
ponent, representing the phonon and photon noises. The very- 
low-frequency drift seen in the real data is added to the TOI. The 
noisy sky TOI is then convolved with the appropriate bolometer 
transfer functions. Another white-noise component, represent¬ 
ing the Johnson noise and readout noise with intensities derived 
in Fig. 26, is also added. Simulated cosmic rays using the mea¬ 
sured glitch rates, amplitudes, and shapes are added to the TOI. 
This TOI is then interpolated to the electronic HFI fast-sampling 
frequency (40 x 181.3737 Hz). It is then converted to analogue- 
to-digital units (ADU) using a simulated nonlinear analogue- 
to-digital converter (ADC). Some 4-K cooler spectral lines are 
added to the TOI. Both effects (ADC and 4-K lines) are derived 
from the measured in-flight behaviour. The TOI finally goes 
through the data compression/decompression algorithm used for 
communications between the Planck satellite and Earth. 

The simulated TOI is then processed in the same way as the 
real mission data for cleaning and systematic error removal, cal¬ 
ibration, destriping, mapmaking, and power spectrum computa¬ 
tion. The systematic effects added to the signal are simulated, 
with the same parameters used by the TOI processing pipeline, 
but no pointing error is included (unlike the simulations of planet 
crossings used for beam and focal plane reconstructions). 

The end-to-end simulations have been used primarily to 
characterize our understanding of the properties of the noise in 
the final maps, such as the level and correlations in the noise 
due to known systematic effects. Figure 27 compares the PSD 
of noise generated by the end-to-end simulation, showing good 
agreement over a wide band: residuals from undetected glitches; 
long time-scale thermal drifts; and filtering of detector and elec¬ 
tronics readout noise. Paper 2 further compares the noise prop¬ 
erties of the end-to-end simulation maps with those of the real 
data. 



10 '^ 10 '^ 10 '^ 10 ° 10 ^ 


Frequency [Hz] 



Figure 26. Two examples of the 4-component noise model fits 
to measured noise PSDs. The underlying noise estimation for 
these fits is described in Planck Collaboration XII (2015). It 
measures the sample auto-covariance of signal-subtracted and 
masked data to overcome effects of gaps and signal residuals. 
The blue curve is the measured PSD averaged over 30 pointing 
periods. The red curve is the 1// component, deconvolved with 
the bolometer transfer function. The black curve is the electronic 
white noise component and the grey curve is the photonic white 
noise component. The magenta curve is the measured ajb corre¬ 
lated component. The green curve is the sum of all noise model 
components. Vertical blue lines indicate the region where the fit 
was made and the vertical black lines indicate the 10 lowest spin 
harmonic frequencies. 


In future releases, the end-to-end simulation will be used to 
characterize the filtering of the sky signal due to the analysis 
pipeline as a function of angular scale. 

6. Conclusions 

This paper has described the major improvements in the HFI 
time-domain processing and beam analysis in the new 2015 
Planck data release. Once the ADC effect is taken into account, 
the bolometer gain stability is at the level of one part in a thou- 
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Figure 27. The PSD of the true noise for bolometer 143-la as 
compared to that generated using the end-to-end simulation. 


sand and the noise is stable within 1 % during the whole mission 
(29 months of sky survey). The '^He-JT 4-K cooler lines have 
been reduced to an insignificant level at the cost of removing 
4 % of the data. The beams are reconstructed out to 100' radius 
at a precision that gives better than 0.1 % errors in the effective 
beam window function. Noise and known systematic errors in 
the beam measurement are negligible in the analysis of the CMB 
data. The characteristics that can be best appreciated at the map 
level are described in Paper 2. Overall, the instrument CMB sen¬ 
sitivity is 13.3 pKcMB Systematic effects have been studied 
individually via dedicated pipelines and globally via end-to-end 
simulations. The main sources of systematic effects are: ADC 
corrections; glitch impacts; 100-mK base-plate temperature fluc¬ 
tuations, 4-K cooler lines, and the bolometer temporal transfer 
function. Their impacts have been measured in terms of trans¬ 
fer functions, noise, and scanning beams. The systematic errors 
here are considered mainly in the context of angular scales with 
{ > 10; as shown in Paper 2, systematic errors remain in the 
data and particular affect the largest scales. A more complete 
description of the instrument and systematic errors in the data 
on the largest scales, f < 10, will be described in a forthcoming 
publication. 

Even with these significant changes in the processing of the 
data and the understanding of the instrument, there is room for 
further improvement. The ADC correction can be improved with 
better modelling of the raw signal and an updated model of 
the ADC code errors. The time response reconstruction method 
which was very effective here for 353 GHz can also be applied 
to 100-217GHz. The cooler lines can be corrected to greater 
precision, allowing more samples of the TOI to be retained for 
mapmaking. All of these corrections in the low-level process¬ 
ing will yield further improvements in the polarization measure¬ 
ments from HFI. 
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Appendix A: HFI TOI and beam product description 

Here we summarize the HFI-specific TOI and beam productes 
that are part of the Planck 2015 data release. A complete de¬ 
scription is given in Planck Collaboration ES (2015). 

All the data products are delivered as FITS files containing 
one or more binary table extensions. Non-data products, such as 
codes, are delivered as tar files. 

A.1. The instrument model 

Formally the “Reduced Instrument Model” or RIMO (since it 
is a subset of the full model maintained internally by the DPC) 
this file contains instrument parameters and some other product- 
related information. 

- Detector parameters: central frequency; beam size, shape, 
position, and orientation in the focal plane; typical noise 
level; calibration factor; and polarization parameters for the 
PSBs. 

- Map-level parameters (for the full-mission, full-channel 
maps only): effective frequency; beam characteristics; and 
noise level. 

- Detector and compound spectral response profiles, for each 
bolometer, for each channel, and for each detector set 
(“Detset”). 

- Detector noise spectra. 

- Beam window functions and the first five error eigenmodes 
for the full channel maps and their cross-products, and for 
the Detset and single (unpolarized) bolometer maps and their 
cross-products. 

- Beam correlation matrices, for the beams listed above. 

Note that only the 50 valid bolometers are included. 

A.2. TOI products 

Clean and calibrated signal TOI is provided in six files, one 
per frequency channel, for each operational day. Each file con¬ 
tains the on-board time (OBT) and associated “global” flags 
(unstable pointing, dark correlation, first and second half-ring), 
and the calibrated signal of all bolometers of that channel (in 
Kcmb for the 100 to 353 GHz channels and MJy sr ' for the 545 
and 857 GHz channels) each with a “local” flag (not-valid data, 
glitch, galactic plane, strong source, etc.). 

The term “clean” means that all known systematic effects, 
including the solar and the orbital dipoles, have been removed, 
and the data, if destriped and projected onto a map using the 
data processing centre (DPC) tools, would yield the DPC maps. 
The flagged regions of the TOI (glitches, big planets, decorre- 
lated dark bolometers) have been replaced by an estimate of the 
signal determined by reading the maps at the pointing position, 
thus they are easily recognized as regions of lower noise. This 
was done to preserve the continuity of the flagged timelines, a 
necessary step for processing them with Fourier techniques. 

In addition to the signal TOI, a pointing TOI is also deliv¬ 
ered. It has the same basic structure as the signal TOI, i.e., one 
file per frequency channel per operational day, but there are three 
columns per detector, corresponding to Galactic coordinates (p, 
9, ip for each sample, i.e., the longitude, the colatitude (= latitude 
-I- 90°), and the roll angle. 

The total volume of these data is around 30 TB, of which 
90 % is the pointing data. 


A.3. ROI products 

Ring-ordered information (ROI) products contain a unique data 
point for each stable pointing period, or ring. In order to repro¬ 
duce the maps delivered by the HFI, some baseline or destriping 
offsets must be applied. These offsets are computed for each de¬ 
tector and for each ring using the full mission data, as described 
in Paper 2. Three sets of offsets are provided: one for the full 
rings, and two for each half-ring. These offsets are given in a 
large table which is Airings rows xAdets columns = 26766 x Ab, in 
which each cell contains a vector with the three offsets, where Nb 
is the number of bolometers with a delivered TOI. These values 
are given in the second extension. The first extension contains 
the indices of the first sample of each ring (and also the ESA 
pointingID and the ring start time) to indicate where each offset 
value has to be applied. 

A.4. Software 

- Colour-correction / unit-conversion code: this is a suite of 

IDF scripts to manipulate the bandpass profiles contained in 

the RIMO to determine colour corrections and unit conver¬ 
sions. 

- Fikelihood code: see Planck Collaboration XI (2015). 

Appendix B: Scanning beam pipeline 

This section describes the special processing applied to planet 
data, leading to the scanning beam maps, emphasizing changes 
from Planck Collaboration VII (2014). We use four observations 
of Saturn and four to five observations of Jupiter for the 100- 
353 GHz channels.^ For these channels, we chose to trade the 
excellent cross-scan coverage of the first Mars observation for 
a higher signal-to-noise ratio. This choice introduced complica¬ 
tions into the pipeline in order to handle data from observations 
separated in time by months, as well as vastly different planet 
brightnesses. 

The submillimetre channels (545 GHz and 857 GHz) are 
driven nonlinear by Jupiter and Saturn, so we use two obser¬ 
vations of Mars to build the beam near the peak of the main 
beam, supplemented by Saturn and Jupiter data to measure the 
near sidelobes. 

The representation of the main beam’s peak and its closest 
sidelobes (to approximately 30') uses B-spline functions (Planck 
Collaboration VII 2014), while the diffraction pattern beyond 
this radius follows an analytic function A9^^, where A is an am¬ 
plitude factor derived from the mirror edge taper (Tauber et al. 
2010 ) and 6 is the angular distance from the main beam axis. 
For beams using Saturn data in the peak, we apply the correction 
for the planet’s finite size at the effective beam window function 
level (Sect. 4.2). Since the scanning beam maps are evaluated 
in the detector frame, the scanning direction is always the same, 
and because the spacecraft rotation rate is nearly constant, tem¬ 
poral residuals can be spatially constrained in those maps. We 
apply a specific treatment for the time-response residual region, 
to improve sensitivity to the residual; specifically, B-splines are 
used instead of the diffraction model. See Sect. 3.4.1 for esti¬ 
mates of transfer-function residuals that cannot be accounted for 
in the scanning beam, due to the finite size of the reconstructed 
maps used for window function estimation. 

^ HFFs dilution cooler began to warm during the fifth observation 
of Jupiter, making only four observations available for the 100 GHz 
bolometers and the 217 GHz SWBs. 
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B.1. Planet data processing 

Planet data from the standard TOI processing pipeline cannot 
be used to recover the scanning beam to sufficient accuracy 
(< 1% at the window function level). We apply a suite of post¬ 
processing operations to the planet data in order to reduce the 
contribution of identified systematic errors to the bias in the 
scanning beam representation and the window function estimate. 
Simulations show that the additional processing applied to plan¬ 
etary data reduces the effective window function bias to less than 
the measurement errors. The following operations are applied 
to the planet data and are specihc to the beam reconstruction 
pipeline: 

- background subtraction; 

- planet proper motion pointing correction; 

- destriping; 

- pointing offset correction per planet crossing; 

- merging of different calibrators (renormalization); 

- additional glitch residuals flagging. 

The background subtraction uses the full-mission, full- 
frequency science maps (the planet crossings are masked). A 
cubic-spline interpolation is performed from the pixel centres 
to the planet crossing pointing (Fig. B.l). 




Ecliptic longitude [°] Ecliptic longitude [°] 

Figure B.l. Hit count map (left panel) of the corresponding 
background map (right panel) of the first Saturn scan, for 
bolometer 143-5, in Kcmb- 


The proper motion of the planets is corrected using 
ephemerides from the JPL Horizons package® (Giorgini et al. 
1996) to rotate the pointing into the Planck spacecraft rest frame. 

The destriping removes long-timescale baseline drifts. Here 
we subtract a baseline offset for each circle (60s timescales) 
rather than for each ring (30-60 min time scales) as in the map¬ 
making. We estimate the baseline offset in a range from 3° to 
6 ° away from the beam peak as the median of the samples prior 
to the planet crossing, to minimize the bias from time-response 
residuals. We then smooth the overall trend per circle using 
a sliding window with a 40-circle width, as seen in Fig. B.2. 
Figure B.3 shows an example of the destriping procedure. 

During the destriping process, we flag and remove the en¬ 
tire circle if one of the following criteria is met in the baseline 
timeline: 

® http://ssd.jpl.nasa.gov/?horizons 


- the number of flagged samples exceeds 25 % of the total 
samples in the circle; 

- there is an event exceeding 3.3 times the estimated local rms 
(residual glitch biasing the estimate of the baseline); 

- the number of positive and negative samples is not balanced, 
such that their relative contribution in the total number of 
used samples differs by more than 20 %. 

An additional 5 % of the planet data are flagged during this pro¬ 
cedure, largely in the main-beam area. 



Circle timeline (arbitrary origin) 

Figure B.2. The baseline for each circle for bolometer 100-3b 
and the hrst scan of Mars (black line), and smoothed with a 40- 
circle sliding window (red line). 


We expect an error of a few arc seconds in the pointing cen¬ 
troid of each planetary observation, given the current pointing 
model (Planck Collaboration I 2015). We apply a re-centreing 
procedure to each planet crossing by fitting a beam template to 
each observation, where the free parameters include the pointing 
shifts in both directions (co-scan and cross-scan). The pointing 
of each planet crossing is corrected with these offsets. The planet 
data processing is an iterative process, with a loop over the beam 
template used. For the hrst iteration, the template is a prior ver¬ 
sion of the beam. 

Each planet crossing has a different peak brightness, so we 
renormalize the planet signals in order to merge them into a com¬ 
mon beam map. We interpolate the beam template to the point¬ 
ing of each crossing to take into account the non-uniform cover¬ 
age of the calibrator observation and perform an azimuth angle 
average (see Fig. B.4) of the samples. We then recover an am¬ 
plitude factor by htting the azimuthally-averaged planet data to 
an azimuthally averaged interpolated template. For Jupiter data, 
only the portion of the prohle not affected by nonlinearity is used 
in the htting procedure. Since the same template is common to 
all calibrators, this procedure, by construction, forces agreement 
between the planet data. Figure B.5 shows the renormalized pro- 
hles. We apply an interpolation on the two-dimensional beam 
map template to take into account the change in the coverage per 
scan, distorting the uniformly distributed azimuthal template, 
since beams are asymmetric and the sampling is sparse during 
a planet observation due to the Planck scan. The spread at low 
radius between different scans and planets is dominated by this 
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Sample timeline (arbitrary origin) 

Figure B.3. TOI for bolometer 100-3b during the first scan of 
Mars, shown both before (upper panel) and after (lower panel) 
subtraction of the baseline (Fig. B.2). 


effect, while the spread at high radius is mainly attributed to the 
different signal-to-noise properties of each calibrator. 

The performance of the standard TOI deglitching is degraded 
by the high signal gradients of planet scans, relative to CMB 
data; with no additional flagging, glitch residuals create a bias of 
several percent in the effective beam window function. We there¬ 
fore add an extra stage of glitch flagging to the planet data pro¬ 
cessing that excludes samples exceeding 3 times the estimated 
local noise rms in the low gradients region of the main beam, 
and 3 times the expected rms of high-frequency pointing errors 
in the strong gradient regions of the main beam (see Fig. B.6). 
The fact that the excluded samples are vertically aligned when 
projected on the cross-scan axis indicates that they are likely to 
be within the same scanning circle, correlating them in time. The 
CO- and cross-scan projection indicates they are localized in the 
beam peak region, where the performance of the standard glitch 
removal process is seriously affected. 

The planetary data processing converges in five iterations to 
within the Monte Carlo error bars (see Fig. B.7), and is not de¬ 
pendent on the template used in the first iteration. 



Distance from beam centroid [arcmin] 


Figure B.4. Renormalization of the beam. The black stars 
show the beam azimuthal profile template used to renormal¬ 
ize the data. The pointing information of the scan is reflected 
through the values of the planet scan’s template, since the two- 
dimensional scanning reference beam map is interpolated to the 
transit pointing before the azimuthal averaging. The renormal¬ 
ization convention is arbitrary, and is chosen so that the template 
resulting from uniform azimuthal coverage is set to the value 1 at 
1' from the beam centroid. The blue curve traces the positive data 
from the planet scan azimuthal average after the renormalization 
procedure; the yellow dashed curve shows the negative data con¬ 
verted to positive values for plotting purposes. For Jupiter, only 
data below the conservative threshold are used in the renormal¬ 
ization procedure (see table B.l) to avoid biases from the finite 
planetary disc size and detector’s nonlinearity response effects. 


<n 



Distance from beam centroid [arcmin] 


Figure B.5. Renormalized planet data (positive values). Planet 
scans are prepared for the merging into a single data timeline us¬ 
ing the azimuthal profile template process. For illustration pur¬ 
poses, we also show Neptune and Uranus scans. 


B.2. Hybrid representation for main beams 

We use Saturn data to measure the main beam and Jupiter to 
measure the closest sidelobes for all channels, except 545 and 
857 GHz where the peak is described with Mars data. The tran¬ 
sition regions consists of a linear weighting scheme between the 
data sets. The transition region is defined as the span of signal 
samples that are lower than 9 times the estimated noise rms of 
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Figure B.6. Additional flagging process in the beam pipeline. Data for samples near the peak response (marked with red crosses in 
the upper panels) are rejected when the residual timeline (planet data minus the template, blue points in the lower panels) is either 
beyond the threshold set by 3 times the rms noise (inner red lines), or beyond the pointing error threshold (outer red lines). 


Table B.l. Jupiter threshold values as a percentage of peak max¬ 
imum amplitude. 


Band [GHz] 

100 

143 

217 

353 

545 

857 

Threshold [%] 

50 

33 

20 

10 

1 

1 


the peak planet (Saturn or Mars) and greater than the Jupiter 
threshold. The Jupiter threshold is set such that we keep only 
data where nonlinearities of the detectors and the finite size of 
Jupiter’s disc have negligible effects compared to the noise un¬ 
certainties. Table B.l summarizes the Jupiter threshold levels in 
peak maximum amplitude. They have been set to conservative 
values. 

We further spatially average the combined Jupiter data in re¬ 
gions where the signal is lower than 3 times the estimated Jupiter 
rms noise. Physical optics simulations show that we should not 
expect high-spatial-frequency variations in the beam in those re¬ 
gions. The regularization filter is a square boxcar average, 0.'5 
on a side for 100-353 GHz, and 10" on a side for 545-857 GHz. 

For the final beam model, we fit the B-spline coefficients 
to the renormalized planet timelines and the filtered Jupiter 
data, using a least squares estimator. As opposed to Planck 
Collaboration VII (2014) (in particular, their equation C.4), no 


additional smoothing criterion is added to the scoring func¬ 
tion in the inversion process, since the signal-to-noise ratio of 
four Saturn crossings is high enough to prevent bias in the B- 
spline functions coming from high-frequency noise. The inver¬ 
sion process uses the Cholesky method from the LAPACK module 
(Anderson et al. 1999). Table B.2 summarizes the control point 
(knot) parameters. The knot locations are tuned so that the fi¬ 
nal bias at the window function level is included in the error 
budget. The geometrical shape of the B-spline functions defines 
the properties of the intrinsic low-pass filter acting on the data, 
while preventing the reconstruction being too sensitive to high 
frequency errors. The parameter choice relies on the data quality 
and the coverage of the combined planet scan. 

In the transition area describing the near sidelobes, we fur¬ 
ther process the B-spline representation in order to capture the 
spatial variations of this intermediate regime while allowing a 
smooth transition into the analytic diffraction pattern. To do this, 
we build azimuthally-averaged profiles in a number of ranges 
of azimuth angle, and propagate these templates with a spline 
interpolation to the final map. The number of azimuthal tem¬ 
plates is constrained by the expected asymmetries and secondary 
lobes derived from optical predictions in this regime of the main 
beam. Table B.3 summarizes the number of angular regions at 
each frequency band. Since the azimuthal templates are equally 
distributed in angle, the number of templates sets the separation 
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Table B.2. B-spline knot spacing for the scanning beam model at each frequency band. 


Band [GHz] 

100 

143 

217 

353 

545 

857 

Knot spacing [arcmin] 

2 

1.5 

1.25 

1.25 

1.25 

0.8 



Figure B.7. Convergence process of the planet data post¬ 
treatment for a 143 GHz channel. The first guess (iteration 0) is a 
beam built using a different (obsolete) recovery of the temporal 
transfer function. The figure illustrates how sensitive the recon¬ 
struction is to any change in the time response and in this case 
how it readjusts through the iterations. The dash-dot lines are the 
standard deviation of the spread from the associated Monte Carlo 
ensembles. The number of iterations was chosen so that the iter¬ 
ation variation within the last loop (iterations 4 and 5) is conser¬ 
vatively within the error bars for all channels from 100 GHz to 
353 GHz. 


between each angularly averaged template, and hence the fidelity 
of the reconstruction to the azimuthal asymmetries of the main 
beam in this regime. 

Figure B.8 shows the azimuthally-averaged profiles as a 
function of the radius from the beam centroid. The differences in 
the templates beyond a 5' radius show the optical asymmetry of 
the beam. In addition, the templates close to 90° azimuth angle, 
shown in blue, show time response residual effects. In this re¬ 
gion, the initial B-spline evaluations are kept with no additional 
azimuthal averaging, in order to better reproduce the uncertain¬ 
ties in the time response of the detector and readout electronics. 
The templates are only used in the regions where Jupiter’s rms 
noise is not negligible, down to the assumed symmetric diffrac¬ 
tion pattern dominant regime. 

The intermediate sidelobes beyond 30' from the beam cen¬ 
troid are dominated by a diffraction pattern and are mod¬ 
elled from optical predictions. The predictions are compati¬ 
ble with Jupiter azimuthally-averaged profiles (see, e.g., Planck 
Collaboration XXXI2014). We make the assumption that the in¬ 
termediate sidelobe behaviour is azimuthally symmetric and we 
keep the B-spline representation in the transfer function residual 
region (temporal causal part of the beam) to properly account for 
temporal transfer function residuals. 



Figure B.8. The 32 radial profiles (averaged over a range of az¬ 
imuthal angles) for the 143-5 beam. The azimuthal angles are 
defined so that 90° corresponds to the profile closest to the resid¬ 
ual time response feature (dark blue). 


Appendix C: Time Response 

C.1. Measurement of 353 GHz bolometer time constants 

As for the other HFI detectors, the 353 GHz bolometer thermal 
time response is described in the Fourier domain as the sum of 
single-pole low-pass filters, each parameterized by a time con¬ 
stant and an amplitude. However, the measurement procedure 
followed a different process than for the 100-217 GHz bolome¬ 
ters (described in Sect. 3.4), jointly solving for the optical beam 
and time response. When applied to 353 GHz bolometers, the 
nominal time response procedure created scanning beams with 
strong (several percent of the peak amplitude) non-optical asym¬ 
metries in the scanning (also causal in time) direction, extend¬ 
ing beyond 30' from the beam centroid. This indicated that the 
recovered time constants from the previous method were far 
from the true ones, distorting the corresponding scanning beams. 
Therefore, a new method was developed. 

The procedure described here, created for the 353 GHz time 
response, will potentially be used in the future at all frequencies, 
due to the improved recovery of the very low frequency portion 
of the time response, as well as better reconstruction of the scan¬ 
ning beam close to the centroid. 

C.1.1. Normalization of the glitch slow component 

The stacked glitch template used for glitch removal (Sect. 3.2 
and Fig. 10; see also Planck Collaboration X 2014) is built using 
data filtered with a 3-point finite impulse response filter.^ Hence, 
the time response inferred from glitches, and in particular the 
relative amplitude of slow response to fast response, is strongly 
biased and must be renormalized. 


’ The data are convolved with the kernel [0.25,0.5,0.25] in the time 
domain. 
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Table B.3. Azimuth angle averaging parameters for linking the closest sidelobes to the diffraction pattern, around 25' from the beam 
centroid. 





Band [GHz] 




100 

143 

217 

353 

545 

857 

Number of azimuthal templates. 

16 

32 

32 

32 

64 

64 

Angular separation. 

22?5 

11?2 

1U2 

1U2 

5?6 

5?6 

Azimuthal averaging upper threshold [Jupiter noise rms]. 

9 

9 

7 

7 

5 

5 

Diffraction pattern dominant regime [Azimuthal averaged Jupiter profile rms] .... 

3.3 

3.3 

3.3 

3.3 

3.3 

3.3 


C.1.2. Starting parameters and dipole correction 

The time constants and associated amplitudes where r > 300 ms 
are set to the glitch values, but renormalized in overall ampli¬ 
tude using planet data with the scanning beam method described 
in Appendix B. The results of the fitting procedure for the time 
response is dependent on the initial guess because of the strong 
correlation in parameters and because planet data (considered 
only within 100' of the centroid) only weakly constrain the very 
slow portion of the time response (r > 1 s). The initial guess 
of the glitch renormalization amplitude is set such that the cre¬ 
ated power compared to the previous release matches the correc¬ 
tion needed to remove half of the shifted dipole effect (Planck 
Collaboration VIII 2015), assuming that the other half of the 
change will be provided by the change of the fast components 
of the time response. 

Because of the convention we adopted for the global normal¬ 
ization of the time response model, any additional response at 
very low frequencies added to match the long timescale response 
of the glitches changes the amplitudes of the existing time re¬ 
sponse. Thus, for all the other parameters, the starting values are 
the previous-release parameters, with amplitudes adjusted to in¬ 
clude the power introduced by the renormalized glitch template 
in the time response. 


C.1.3. Beam-forming while recovering the parameters 

The simultaneous fitting of an optical beam shape and time- 
response parameters for the 353 GHz detectors follows the pro¬ 
cedure outlined in Planck Collaboration VII (2014), apart from 
slight differences in the distance from the planet centroid where 
the fitting is performed. A symmetrized beam profile along the 
scanning direction is derived from initial-guess time response 
parameters, then updated after each iteration. 

The fit is performed as a loop over two steps, alternating fast 
(10ms < T < 100ms) and slow (100ms < t < 300ms) param¬ 
eter adjustments. For fast parameters, the recovery is performed 
over the range 3' to 70' from the beam peak centroid, while for 
slow parameters the range starts at 20' and extends up to 300'. 
In the first step, the fast parameters are fitted to combined Saturn 
and Jupiter data, with other parameters fixed (intermediate ones 
and glitch template amplitudes). The second step adjusts inter¬ 
mediate parameters and the glitch template amplitude to match 
Jupiter data with other parameters fixed (fast parameters). We 
stop the loop when both of the following requirements are met: 
(1) the symmetrized beam template is stable within 1 % accu¬ 
racy in the peak domain; and (2) the relative reduced is stable 
within 10 in each iteration. Figure C.l illustrates the overall 
process for the 353 GHz channels. The final parameters are listed 


Saturn data 

Jupiter data 

(4 transits) 

j (5 transits) j 


1 1 * 


Data merging 

using the same planet data 
post-processing as in the 
scanning beam pipeiine 


Rectangular data selection : 
[-0.5', 0.5'] cross scan width 
[-30', 300'] in scan width 
with origin (0,0) set on the pianet peak 


Coliapse into a single in-scan 
oversampled timeiine 
using B-Spline representation 


'Opticai' beam model 
(guess of the underlying opticai beai 
free from time constants effects, 
in scan symmetrized) 



Beam timeline 

convolved with 
the time response 


initial parameters 

Beam model; previous release + close side lobes! 
Temporal transfer function : previous reiease h 
re-normaiized glitches tempiates + 
dipoie shifts constraints 


Fast parameters recovery 
using planet data : 

[2', 70'] from beam centroid 

Slow parameters recovery 
using pianet data : 

[20',300'] from beam centroid 


Figure C.l. Schematic diagram of the time-constant fitting pro¬ 
cess for the 353 GHz channels. The main difference from the 
CMB channels is the fact that the longest time constants from 
the glitch templates that cannot be constrained by the planet sig¬ 
nal are renormalized so that their relative contribution is injected 
properly in the temporal model, and the fast part of the time re¬ 
sponse is re-evaluated to reduce the co-scan distortion and in¬ 
duce a beam closer to optical. This is particularly significant 
for PSB pairs, where incorrectly deconvolved fast time constants 
can artificially induce beam mismatch and therefore spurious po¬ 
larization. 


in Table C.l. Notably, this solution requires the inclusion of two 
extra time-constants and associated amplitudes. 


^ APC, AstroPaiticule et Cosmologie, Universite Paris Diderot, 
CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris 
Cite, 10, rue Alice Domon et Leonie Duquet, 75205 Paris Cedex 
13, France 

^ African Institute for Mathematical Sciences, 6-8 Melrose Road, 
Muizenberg, Cape Town, South Africa 

^ Agenzia Spaziale Italiana Science Data Center, Via del Politecnico 
snc, 00133, Roma, Italy 

'* Aix Marseille Universite, CNRS, LAM (Laboratoire 
d’Astrophysique de Marseille) UMR 7326, 13388, Marseille, 
France 

^ Astrophysics Group, Cavendish Laboratory, University of 
Cambridge, J J Thomson Avenue, Cambridge CB3 OHE, U.K. 

® Astrophysics & Cosmology Research Unit, School of Mathematics, 
Statistics & Computer Science, University of KwaZulu-Natal, 
Westville Campus, Private Bag X54001, Durban 4000, South 
Africa 
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Table C.l. Parameters for time response models that are deconvolved from the data. 


Bolometer 

a\ 

[10-3] 

Tl 

[ms] 

cii 

[10-3] 

Tl 

[ms] 

as 

[10-3] 

U 

[ms] 

a4 

[10-3] 

U 

[ms] 

as 

[10-^] 

Ts 

[ms] 

06 

[10-'] 

n 

[ms] 

d-j 

[10-'] 

Tl 

[ms] 

^stray 

[ms] 

^ phase 

[ms] 

100-la 


394 

10.0 

537 

20.9 

54.6 

52.3 

11.6 

246 

2.09 

1375 





1.593 

1.386 

100-lb 


478 

10.3 

457 

19.2 

45.9 

56.3 

13.4 

305 

5.74 

2187 





1.489 

1.386 

100-2a 


483 

6.84 

429 

13.6 

72.8 

43.1 

13.3 

298 

1.28 

1910 





1.321 

1.247 

100-2b 


132 

5.84 

749 

15.1 

104 

43.4 

13.4 

409 

2.26 

2925 





1.378 

1.247 

100-3a 


741 

5.39 

222 

14.7 

29.4 

47.7 

7.88 

327 

0.09 

2113 





1.421 

1.247 

100-3b 


601 

5.48 

348 

15.5 

38.8 

43.5 

9.83 

249 

2.37 

1058 





1.663 

1.247 

100-4a 


410 

8.20 

513 

17.8 

60.6 

54.4 

11.8 

325 

4.41 

1839 





1.251 

1.247 

100-4b 


690 

11.3 

283 

24.3 

20.3 

90.0 

2.26 

295 

4.57 

1924 





1.381 

1.386 

143-la 


816 

4.47 

143 

12.0 

27.7 

39.4 

11.7 

281 

1.03 

1250 





1.418 

1.247 

143-lb 


498 

4.72 

338 

15.6 

129 

55.7 

32.1 

352 

3.44 

1656 





1.485 

1.247 

143-2a 


910 

4.70 

76.4 

17.0 

8.47 

100 

4.36 

290 

1.19 

1370 





1.481 

1.247 

143-2b 


925 

5.24 

51.6 

16.7 

15.2 

49.0 

7.37 

279 

0.63 

1269 





1.461 

1.247 

143-3a 


687 

4.19 

276 

9.56 

27.3 

43.6 

8.34 

308 

1.40 

1499 





1.451 

1.247 

143-3b 


820 

4.47 

131 

13.2 

36.9 

39.9 

10.0 

324 

1.65 

1735 





1.609 

0.832 

143-4a 


911 

5.69 

71.7 

18.9 

10.6 

49.4 

5.38 

318 

1.78 

1786 





1.592 

1.247 

143-4b 


429 

6.06 

508 

6.06 

55.6 

26.7 

7.05 

363 

0.46 

2208 





1.820 

1.247 

143-5 


523 

6.64 

423 

6.64 

38.9 

42.5 

11.4 

322 

3.76 

1425 





2.024 

1.386 

143-6 


535 

5.51 

423 

5.51 

31.5 

38.7 

9.39 

338 

1.53 

2227 





1.529 

1.109 

143-7 


415 

5.43 

565 

5.43 

14.8 

43.3 

4.40 

366 

0.80 

2708 





1.860 

1.386 

217-1 


13.6 

3.46 

955 

3.46 

27.5 

25.0 

2.95 

384 

0.96 

3284 





1.590 

1.109 

217-2 


979 

3.52 

14.0 

26.1 

2.91 

41.7 

2.73 

362 

1.28 

2748 





1.602 

1.247 

217-3 


934 

3.55 

33.7 

3.55 

27.8 

31.7 

4.20 

321 

0.38 

1754 





1.742 

1.247 

217-4 


658 

1.35 

320 

5.55 

17.9 

27.9 

2.74 

402 

0.80 

2920 





1.710 

1.109 

217-5a 


905 

6.69 

79.7 

21.6 

7.68 

66.2 

6.34 

252 

1.79 

1160 





1.573 

1.109 

217-5b 


924 

5.76 

60.9 

18.0 

7.49 

65.8 

6.19 

343 

1.42 

2742 





1.867 

1.247 

217-6a 


872 

6.45 

69.7 

19.7 

46.8 

44.8 

11.3 

270 

0.46 

1226 





1.545 

1.247 

217-6b 


285 

6.23 

667 

6.23 

38.9 

26.8 

8.65 

267 

0.22 

1266 





1.455 

1.109 

217-7a 


344 

5.48 

576 

5.48 

71.6 

25.1 

6.82 

282 

1.77 

1279 





1.515 

1.386 

217-7b 


847 

5.07 

127 

14.4 

17.3 

49.9 

6.19 

348 

2.36 

1787 





1.505 

1.386 

217-8a 


498 

7.22 

441 

7.22 

50.7 

30.2 

9.05 

266 

1.04 

1260 





1.789 

1.109 

217-8b 


512 

7.03 

411 

7.03 

63.5 

27.7 

9.65 

312 

1.22 

2212 





1.731 

1.247 

353-1 


740 

1.23 

187 

6.07 

59.3 

16.1 

106 

71.9 

1.92 

1378 

0.57 

1725 

0.90 

4034 

1.211 

0.911 

353-2 


944 

5.32 

44.9 

19.8 

6.63 

103 



0.86 

879 

2.99 

2611 

0.17 

7084 

3.445 

0.970 

353-7 


686 

0.55 

283 

3.60 

24.5 

14.6 

3.98 

99.9 

2.69 

682 

0.05 

1180 

0.22 

3498 

1.089 

1.247 

353-8 


716 

1.97 

261 

5.95 

17.3 

31.7 



4.25 

169 

0.59 

1052 

1.51 

2943 

1.987 

1.109 

353-3a 


22.7 

3.03 

869 

6.90 

81.2 

22.5 

16.2 

80.5 

10.3 

490 

0.11 

1310 

0.17 

3611 

1.782 

1.247 

353-3b 


218 

2.67 

667 

6.95 

88.9 

19.5 

19.4 

70.0 

5.73 

597 

0.08 

1402 

0.19 

3770 

1.574 

1.109 

353-4a 


463 

1.89 

464 

6.26 

55.5 

22.1 



10.8 

119 

4.94 

1507 

1.69 

3311 

1.479 

1.247 

353-4b 


809 

4.50 

89.0 

15.5 

19.9 

28.9 

73.5 

96.2 

8.68 

326 

0.24 

2181 

0.08 

4611 

1.726 

1.109 

353-5a 


778 

5.94 

160 

12.4 

41.3 

32.4 

14.2 

104 

0.21 

559 

6.59 

1334 

0.08 

3356 

1.571 

1.109 

353-5b 


782 

5.89 

115 

10.8 

27.3 

39.9 

67.3 

9.57 

6.15 

272 

0.65 

1429 

1.51 

3510 

1.893 

1.109 

353-6a 


93.7 

1.27 

834 

5.99 

58.0 

24.2 

11.1 

112 

4.06 

701 

0.20 

1565 

0.15 

3666 

1.998 

1.247 

353-6b 


53.7 

1.24 

911 

5.54 

23.1 

28.8 

6.29 

118 

5.84 

763 

0.07 

1201 

0.18 

3444 

2.895 

1.109 

545-1 


991 

2.93 

7.43 

26 

1.39 

2600 









2.160 

1.109 

545-2 


985 

2.77 

12.8 

24 

2.46 

2800 









1.868 

0.970 

545-4 


972 

3.00 

27.7 

25 

0.78 

2500 









2.222 

1.109 

857-1 


974 

3.38 

22.9 

25.0 

3.49 

2200 









1.765 

1.109 

857-2 


840 

1.48 

158 

6.56 

2.49 

3200 









2.202 

1.247 

857-3 


360 

0.04 

627 

2.40 

11.1 

17 

2.00 

1900 







1.524 

1.263 

857-4 


278 

0.40 

719 

3.92 

1.62 

90 

1.52 

800 







1.490 

0.558 
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